PDE System Approach to Large Extra Dimensions 



Edward Anderson and Reza Tavakol 



Astronomy Unit, School of Mathematical Sciences, 
Queen Mary, Mile End Road, London El 4^S, U.K. 

Abstract 

We explore some foundational issues regarding the splitting of ZJ-dimensional Einstein Field Equations 
(EFE's) with respect to timelike and spacelike {D - l)-dimensional hypersurfaces, first without and then with 
thin matter sheets such as branes. We begin to implement methodology, that is well-established for the GR 
Cauchy and initial value problems (CP and IVP), in the new field of GR-based braneworlds, identifying and 
comparing many different choices of procedure. 

We abridge fragmentary parts of the literature of embeddings, putting the Campbell-Magaard theorem into 
context. We recollect and refine arguments why York and not elimination methods are used for the GR IVP. We 
compile a list of numerous mathematical and physical impasses to using timelike splits, whereas spacelike splits 
are known to be well-behaved. We however pursue both options to make contact with the current braneworld 
literature which is almost entirely based on timelike splits. 

In our study of timelike splits, we look at the Shiromizu-Maeda-Sasaki braneworld by means of reformula- 
tions which emphasize different aspects from the original formulation. We show that what remains of the York 
method in the timelike case generalizes heuristic bulk construction schemes. We formulate timelike (brane) 
versions of the thin sandwich conjecture. We discuss whether it is plausible to remove singularities by timelike 
embeddings. We point out how the braneworld geodesic postulates lead to futher difhculties with the notion 
of singularities than in GR where these postulates are simpler. 

Having argued for the use of the spacelike split, we study how to progress to the construction of more general 
data sets for spaces partially bounded by branes. Boundary conditions are found and algorithms provided. 
Working with (finitely) thick branes would appear to facilitate such a study. 
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1 Introduction 



A great deal of effort has recently gone into the study of higher-dimensional models as unified theories of gravity 
and fundamental matter fields. The original interest in such models in Kaluza-Klein (KK) theory ^ was revived 
by the important role of dimension 11 in supergravity |2] and more recently by superstring theory and M- 
theory which favour spacetimes of dimension 10 and 11 respectively. The usual (KK) argument for the 
apparent 4-dimensionality of spacetime is that the extra dimensions are compactified i.e curled up sufficiently 
small so as not to conflict with observation. However, recently much attention has been placed instead on string- 
theory-inspired models with large extra dimensions [HI El 13 13 El ^1 in which most of the physics is confined 
or closely-bound to a lower-dimensional braneworld surrounded by a higher-dimensional bulk. 

Though motivated by superstring and M-theories, many of these models El d IH1 are in fact formulated 
within the framework of (higher-dimensional) general relativity (GR). Among these are: 1) the second Randall- 
Sundrum scenario [7] in which the graviton is tightly-bound to the brane by the curvature due to the warping 
of the bulk metric, which is pure anti-deSitter (AdS). 2) A more general scheme of Shiromizu, Maeda and 
Sasaki (SMS) in which the 4-dimensional Einstein field equations (EFE's) are replaced by 4-d "braneworld 
EFE's" (BEFE's), which are not closed since there is a 'dark energy' Weyl tensor term, knowledge of which 
requires solving also for the bulk. The interpretational difficulties due to this constitute the 'Weyl problem'. 
The other main distinctive feature of SMS's BEFE's is the presence of a term quadratic in the braneworld 
energy-momentum, which arises from the junction condition used jl2| IHj. 

Two important questions arise in consideration of such models. Ql: how should such models be built 
and interpreted consistently within the framework of GR? This would require a careful underlying choice of 
conceptually-clear general framework, in the sense we discuss below. Q2: what exactly is the connection between 
such models and the underlying string/M theory? More precisely, to what extent can the agreements or otherwise 
of predictions of such models with observations be taken as support or disagreement with such theories? Here we 
concentrate on the first question and make a comparative study of the general schemes that have been employed 
in the literature in order to construct the bulks which surround branes. 

Such a comparative study requires a sufficiently general common language. We use that of p.d.e problems, 
which consist of both the p.d.e system to be solved in some region of a manifold, and data prescribed on (portions 
of) the boundary of this region. For (r, s) a spacetime of dimension n = r + s with s time dimensions, we denote 
the problem involving the addition of an extra dimension by (r, s; e), where e = 1 if the new dimension is 
spacelike or e = —1 if it is timelike. This is the generalization of the GR Cauchy problem (CP) jl3l 1141 [T^ I16j 
based on the arbitrary (r, s; e) generalization of the Arnowitt-Deser-Misner ^2] split of the metric (Sec 2.1) 
and hypersurface geometry (Sec 2.2). Its simple signature-independent features are pointed out in Sec 2.3 and 
the crucial dependence of many of the harder features on the usual CP signatures s = 0, e = — lis discussed in 
Sec 2.4. As is well-known, the GR CP presupposes the existence of the data. Thus one is in fact considering a 
two-step process, the other step being the construction of the data on the (r, s) manifold i.e the generalization 
of the GR initial value problem (IVP) [IHlinil^- This is discussed in Sec 3. 

Our framework permits a profitable look at a number of recent topics. The aim is to consider GR-based 
models containing thin matter sheets such as branes or domain walls. We shall compare two broad schemes that 
have been proposed to construct bulks: the (3, 1; 1) construction [211 1221 OH] starting from information on a 
(3, 1) spacetime hypersurface (usually the brane) and the (4, 0; -1) construction U4, 25 starting by the construc- 
tion of data on a (4, 0) spatial hypersurface. However, first we emphasize that one should grasp the fundamental 
arguments and results which before the specializing to the thin matter sheet models. A few ideas about the 
general unspecialized framework recently arose in the literature on (generalizations of) the Campbell-Magaard 
arbitrary-signature embedding theorem |26 | \27 \ 1 ^ I29j (see Sec 2.3). However, we find it far more profitable 
instead to adopt the generalized GR CP-IVP point of view since this literature is by far more developed and 
thus a far greater source of well-documented pitfalls and carefully thought-out techniques which avoid them. 
We identify the embedding step of the Campbell-Magaard theorem with the well-known signature-independent 
parts of the GR CP in Sec 2.3. But the harder s = 0, e = — 1 specific parts of the GR CP strongly suggest that 
(4, 0; -1) schemes should be favoured over (3, 1; 1) ones on very general grounds: well-posedness and causality 
(Sec 2.4). Furthermore, Magaard's data construction method ^\ (Sec 3.1) does not compare favourably with 
York's data construction (Sec 3.2), and its application to s = 1, e = 1 has further conceptual difficulties. In 
this light we look at the extent to which York's method is adaptable to s = 1, e = 1, and also consider the thin 
sandwich method |3I3I^ iii this context (Sec 3.3). 

We then introduce thin matter sheets in Sec 4, and study the (3, 1; 1) schemes with thin matter sheets, 
recollecting the derivation of the junction conditions (Sec 4.1), and showing how the SMS formulation (Sec 4.2) 



may be reformulated in a large number of ways using geometrical identities that interchange which terms are 
present in the braneworld EFE's. This is illustrated by the formulation in Sec 4.3 which directly parallels the 
GR CP formulation and thus makes no explicit use of the Weyl term, and by formulations in Sec 4.4 in which the 
quadratic term has been re-expressed entirely in terms of derivatives off the brane. These formulations are used 
to clarify a number of aspects of the 'Weyl problem', in particular to argue for the study of the full brane-bulk 
system. Further aspects of such formulations are discussed in Sec 4.5. The implications of Sec 2.4 for thin matter 
sheets, the heuristic use and limitations of the (3, 1) York method with thin matter sheets, and formulations of 
thin matter sheet thin sandwich conjectures are discussed in Sec 4.6. 

Having gathered arguments against using (3, 1; 1) schemes in Sees 2.4 and 3, we further consolidate this 
viewpoint by arguing in Sec 5 against a suggested virtue of (3, 1; 1) schemes [211 IS]: their use to remove 
singularities. In addition to using the framework of Sec 2, the corresponding study of geodesies is also used here, 
which presents further conceptual complications relevant to thin matter sheet models. These follow from the 
(3, 1) geodesies not usually being among the (4, 1) geodesies. 

In Sec 6, we continue the study of the (n, 0; -1) scheme favoured by our arguments, in the presence of both 
thin and thick (i.e finitely thin) matter sheets. We begin in Sec 6.1 by providing a hierarchy of very difficult 
general thin matter sheet problems in which the main difficulties stem from low differentiability and details about 
the asymptotics. Within this class of problems we identify how the far more specific scenarios currently studied 
emerge as more tractable cases, and thus identify which as-yet unjustified assumptions such studies entail. The 
IVP step, of use in the study of how braneworld black holes extend into the bulk (the 'pancake' versus 'cigar' 
debate [25)' involves fewer of such assumptions. Thus in this paper we restrict attention to the (n, 0) data 
construction problem (Sec 6.2), which we apply to thin matter sheets in Sec 6.3 and more straightforwardly to 
thick matter sheets in Sec 6.4.. Sec 7 discusses Q2 and contains the conclusions of this paper as regards Ql. 

2 Embeddings and the GR Cauchy Problem 

In this section, we give splittings of EFE's with respect to a hypersurface of codimension 1. This results in 
2 systems of equations that are treated as constituting two separate steps. The first, embedding or evolution 
step, will be dealt with in this section and the second, data construction step, will be presented in Sec 3. The 
embedding step includes the GR Cauchy problem (CP) for s = 0, e = —1. We work with general (r, s; e) since 
many applications of interest have s = 1 and e = 1. Our approach is to treat the general case as far as possible, 
while taking careful note of the many important roles played by the signatures e and especially s (whereas a 
change in the dimension n plays only a mild role). 

An important point regarding the well-known analytic p.d.e system approach employed here is the identifi- 
cation of the input to be prescribed and the output to expect. With the advent of schemes to construct bulks 
[^ [^ [2H1 [m [33 [SSj ! the old questions^ regarding these prescriptions have become relevant in a new context. 

2.1 Splitting up the Metric 

We begin by recalling the arbitrary (r, s; e) Arnowitt-Deser-Misner (ADM) jl7j type split w.r.t the (r, s) 'initial' 
hypersurface Tq of the higher-dimensional metric gcD^ along with the split of its inverse g'~'^ , 



The geometric interpretation of (and further notation for) the (n, 0, -1) and (r, 1, 1) splits discussed in this 
paper are exhibited in Fig. 1. For the specific ADM case, the added dimension (of signature e = —1) is time, 
and the lower-dimensional hypersurface is 3-dimensional and entirely spatial: (r, s) = (3, 0). The EFE's are 
then decomposed into 4 constraints and 6 evolution equations (Sec 2.3). The GR IVP and CP are defined w.r.t 
this split. The GR IVP is the solution of the 4 constraints [on some region Sq of some (3, 0) initial hypersurface 
So] for initial data (hat, Kab), where Kab is the extrinsic curvature 



^Following from the original (3, 0; -1) work on GR by Lichnerowicz |36l and (Choquet)-Bruliat |37l 1381 113| (see Sec 2.4 and 
3.1-2), Wheeler |39| asked questions about prescriptions, which stimulated the York and thin-sandwich work (in Sees 3.2 and 3.3 
respectively) . 

^In this paper we use capitals for the higher-dimensional indices involved in a split and lower-case for the lower-dimensional 
indices, qab, hab', Vc, D^; TZabcd, Rabcd', Gab, Gab', Tcd, Td are respectively the higher and lower dimensional metrics; covariant 
derivatives; Riemann tensors (and similarly for Ricci tensors and Ricci scalars); Einstein tensors; energy-momentum tensors. 
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which is a measure of how much the hypersurface is bent in the surrounding spacetime.^ This data, together 
with the 6 evolution equations (for Kab) and the 6 equations for hab which follow from constitute the GR 
CP. This is a (3, 0; ~1) problem for the construction of the spacetime region in which So is embedded. 
It is useful to juxtapose here another metric split, the arbitrary (r, s; e) KK type split 

so as to discuss some of the literature 00]. Note first that for KK theory itself, this split is (3, 1; 1), with 
hcdi CLc, $ held to be independent of the added dimension's coordinate z. This so-called cylindricity condition 
implicitly permits Oc to be interpreted as the classical electromagnetic potential. More recent generalizations 
of this scheme, called 'noncompact KK (NKK) theory' [30], have involved large extra dimensions in place of 
cylindricity, by permitting hcd, clc to depend on z. In some simple instances, Oc is held to be 0. But clearly for 
flc = and Pc = 0, the two splits ^ and Q are identical.'^ So there is already a vast literature on embeddings 
with a large extra dimension since this is what is usually treated in the GR CP. 

In Sec 4, we reinterpret (3, 1; 1) methods such as that of SMS as z-dynamics i.e "dynamics" with the 
coordinate of the new dimension, z, as "dynamical variable".^ Note that such a "dynamical" interpretation was 
of no avail in KK theory proper, for there hcd was taken to be independent of z. Note also that although a way 
of relating SMS's BEFE's and NKK has been pointed out 141 j, we emphasize that it is far more useful as here to 
relate the steps leading to SMS's BEFE's to the GR CP and IVP literature, since this is by far more developed 
and thus a good source of ideas and theorems, on which the critical Sec 2.4 is based. 
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2.2 A Full Account of Hypersurface Geometry 

Gauss' outstanding theorem {R = 2K^ for a 2-sphere embedded in flat 3-space) already provides a good intuitive 
feel both of why matter squared terms arise in brane cosmology and of why intrinsic curvature singularities can in 
principle be removed upon applying an embedding by being balanced out by extrinsic-curvature-squared terms. 
We work with the suitable generalization of this theorem: the (n + l)-dimensional geometry of projections 
w.r.t a (r, s) hypersurface. We use 0,„a... = 0,,,A...ha for projections onto the (r, s) hypersurface T, and 
0...±... = 0,„A...n^ for projections onto the normal. The three projections of the Riemann tensor |42| are 
respectively the Gauss, the Codazzi and the Ricci equations: 

T^abcd = Rabcd - '2eK^[cKd]b^ (4) 



T^±abc = -'^(^D[cKb]a, (5) 

5^1 



T^±a±b = -{SgKab - eDbDaa) + Ka^Kcb- (6) 

a 



The Ricci equation is particularly important in Sec 4.3-4. Contracting these equations gives 

Tlbd - en^b±d = Rbd - e{KKbd - Kb'Kcd). (7) 

Ga± = 7^a± = -e{DbK\ - DaK), (8) 
53K - eD^a 

n±± = KoK. (9) 

a 

The Gauss equation may be contracted a second time to give (using K o K = KijK^^) 

2g_i_i = 2n±±-en = -eR + K'^-KoK. (10) 
The Ricci equation can be used in the contracted Gauss equation to obtain 



-{SgKab - eDbDaa) - KK^b + 2Ka'Kbc + eRab = eUab- (11) 



^We use £(} to denote the Lie derivative w.r.t and a dot to denote the derivative in the direction perpendicular to the 
hypersurface. The hypersurface derivative is given hy50=-^ — £(} where /3 is the 4- vector [a, Pi]. 

*This requires the identification a More generally, a KK split is the inverse of an ADM split with the identifications a <-> ^, 

j3i « — tti. 

^We use \ = z when the new dimension is spacelike and \ = t when it is timelike. We use one more set of apostrophes than is 
necessary to codify the nonstandard- or general-signature use of concepts that are usually considered in the context of a particular 
signature: dynamics, momentum, evolution, the Cauchy problem. 



The trace of this [a hnear combination of H1U|1 and @] is: 



^ {50K - eD^a) -K^ + eR = en- n±± . (12) 



a 

Another linear combination of H10|) and @ is 



a 2 2 

From this and it fohows that 



^ - e^- (13) 



1 

a 



KoK + K"^ 



e{60Kab - habS^K) - DbDaa + habD^a + e [ 2Ka''Kbc - KKab + " ^\ ' hab ]+Gab = Qab (14) 



2.3 Simple Signature-Independent Results for Embeddings 

It is well-known that the general (r, s; e) ADM-type split of the EFE's Gab = Tab gives rise via ()1U[) and ^ to 
the (n + 1) Gauss and Codazzi constraints^ 

-KoK -eR = 2a±± = 2T±± = 2p, (15) 

- eiDbK\ - DaK) = = Ta± = ja, (16) 
together with n(n + l)/2 "evolution equations" w.r.t the "dynamical variable" A 

-{SgKab - eDbDaO) - KKab + 2Ka^Kbc + eRab = elZab = eiSab rrhab ) ^/iafe, (17) 

a P \ n — 1 J n — 1 

obtained via (fTH) . Here p, ja and Sab = Tab (N.B S ^ T) are general matter terms which are usually prescribed 
as functions of matter fields that are governed by usually-separate field equations. A useful equation is the trace 
of (HH) 

6,K = a(K^-eR- !lP±ll) + = a\K o K + (^ " 2)P " ^'^ 1 + ,^2^^ (^8) 

P \ n — 1 J I n — 1 j 

the second equality of which follows from the Gauss constraint H15() . 

Also applying to Kab the decomposition of symmetric second-rank tensors Aab into their trace A and tracefree 
part, = Aab — ^hab, the equations (|T^ and ((TB|) and H18() take on the useful forms 

K^ oK^ - + + 2p = 0, (19) 

n 

DbK^\ - "^—^DaK + eja = 0, (20) 
n 

^{6,K - eD'a) - ^ = K^ o K^ + ^^-^)P-'S . (21) 

We now consider the simpler results of the (3 ,0; -1) GR CP, most of which hold irrespective of both e and 
s. First, there is the two-step nature of the problem for the whole EFE system. 

"Evolution" or embedding step: given as data {hab, Kab) together with coordinate functions p and ja, all 
on some portion Uq of Tq, such that the constraints are satisfied, one is to solve the evolution equations ((21) and 
(|17j) as a prescription for evaluating {hab, Kab) on some portion Ux of a nearby Tx i.e for a small increment of 
the dynamical variable A. 

"initial value" or data construction step: to construct such data by solving the constraints (see Sec 3). 

Second, provided that the constraints hold on some "initial" Uq, they are guaranteed to hold throughout 
the embedding spacetime region. This is most quickly seen by virtue of energy-momentum conservation and the 
Bianchi identities [171, making use of the Einstein tensor projection form of the constraints. 

Third, constraints come hand-in-hand with gauge freedoms. Thus n -|- 1 of the metric components are freely- 
specifiable. For many applications a choice of gauge is required. To begin with we mostly use the normal 
coordinate gauge {a = 1, Pi = 0), and the harmonic gauge |13| I16j (coordinates x"^ such that V^x^ = 0), 

®In ADM's Hamiltonian language for (r, 0; -1) these are known as the Hamihonian and momentum constraints respectively. 



though we use more elaborate gauge choices later. We stress that careful consideration of gauge choice can be 
an important ingredient toward producing prescriptions that are generally robust and physically clear |43j . 

Fourth, by choice of a suitable gauge it can be ascertained by the Cauchy-Kovalevskaya (CK) theorem^ [33] 
for admissible analytic data that there locally exists a unique solution to the "evolution equations" H17|). 

For this and further applications in the GR CP, the harmonic gauge is usually chosen. In the (n, 0; -1) GR 
CP, the use instead of the normal gauge is discouraged because in practice it typically breaks down quickly. To 
see this recall that in normal gauge equation (|21|) may be written as the normal Raychaudhuri equation 

k = K^ oK^ + n±± > 0, 

n 

where the inequality follows from the definition of the strong energy condition (SEC): 

The general Raychaudhuri equation term TZab^^u^ > V unit timelike u^, 
by picking = = t^. To make ()23p manifestly an energy condition, one uses the EFE's to obtain 

Tabu > V unit timelike 

n — 1 

Integrating the differential inequality in ()22() shows that if Kq > 0, \K\ — > oo within a finite range of the 
parameter vr = along t^. This blowup is by definition a caustic and signifies a breakdown of the normal 

coordinates. However, for the (r, 1; 1) embedding n"^ = is spacelike. There is then no good reason for 7^_l_l 
to be always positive, so caustics may be less likely to form. The issue of the longevity of normal coordinates is 
of interest because they have hitherto been used for some of the applications below. Other applications below 
fix /3j = alone. 

An embedding result we discuss below is the Campbell-Magaard (CM) theorem, according to which any 
n-space can be surrounded by n + 1 dimensional vacuum space, where the spaces in question are analytic. This 
is closely related to the above split of the EFE's; in particular it has an embedding step which we consider in this 
section and a data construction step, 'Magaard's method', which we consider in Sec 3.1. The embedding step of 
the CM theorem is simply the result above based on the CK theorem, which is well-known for the GR CP, and 
follows for the other (r, s; e) embedding cases from the well-known insensitivity of the CK theorem to (r, s; e). 
Strictly, the CM result is the vacuum case (Tab = 0) but the proof may easily be extended to any functional 
form of Tab [201 since the use throughout of the CK theorem is totally insensitive to changes in such subleading 
order terms (contributions to Fa). There have been suggestions that such [33] or similar (3, 1; 1) methods 
could be used for the construction of bulks. However, even though the above results in the (n, 0; 1) case of the 
GR CP are protected and extended by a host of further theorems, they do not carry over to the (3, 1; 1) case, 
as we argue in the next section. There are also problems with (3, 1) data construction, as we shall see in Sec 3. 

2.4 The Breakdown of the Cauchy Problem — Embedding Analogy 

In this section we argue by further comparison with the GR CP that the popular practise of building bulks by 
z-dynamics is in general unsatisfactory. For, in the GR CP a large number of important results j37| 146 1 138| I13| 
in. ,47. l5l 1451 148j for the GR CP turn out to be entirely dependent on the lower-dimensional signature s = 0. In 
other words, the choice of methods which properly respect the difference between space and time is absolutely 
crucial. 

To have a sensible problem in mathematical physics, one requires not just local existence and uniqueness of 
solutions but rather well-posedness. Whereas well-posedness always includes local existence and uniqueness, it 
also always includes continuous dependence of the data, without which an arbitrarily small change in the data 
could cause an arbitrarily large immediate* change in the evolution [33j i.e there is no guarantee of physical 

'^Cauchy-Kovalevskaya theorem: Suppose one has A unknowns which are functions of the "dynamical variable" A and of 7 
other independent variables x^. Suppose one is prescribed on some domain M the A p.d.e's of order k of form ^g^^f^ = for the 
_Fa functions of A, x-,, ^'a and of derivatives of 'I' a up to (fc - l)th order w.r.t A, where furthermore in M the -Fa are analytic functions 
of their arguments. Suppose one is further prescribed analytic data ^a(0, x-y) = '"^d^ (x-y), ... , ^3jfc~T^(0' ^7) ~ '"''~^^dA{xj) on 
some piece ?7 of a nowhere- characteristic surface. Then this problem locally possesses precisely one analytic solution. Note 1) this is 
the most basic theorem for p.d.e's. 2) It depends entirely on the functions being analytic, which seriously limits its applicability. 3) 
It is a very general theorem in the sense that it holds for all types of p.d.e's. In this paper, this is reflected by the signature (principal 
symbol)-independence of the applications and by the irrelevance of what analytic function Fa is. 

*We emphasize that the continuous dependence in question (see e.g page 229 of |44|') is an immediate and discontinuous change 
in behaviour, not an issue of chaos or unwanted growing modes (though well-posedness often also bounds the growth of such modes 



(22) 



(23) 



(24) 



predictability from such a problem. For a hyperbolic-type system one further requires a domain of dependence 
(DOD) property to enforce a sensible notion of causality. If the data is known only on a closed achronal set S 
(that is, a piece of a spacelike hypersurface) , then the evolution can only be predicted within a region I?~*"(S') 
(the future DOD), defined as the set of points such that all past-inextendible causal curves through each point 
[represented by the 7's leading to the point p in Fig. 4a)] intersect S. In our view causality can effectively be 
studied only in settings where the dynamical variable is time. One reason for this is that given information on 
an arbitrarily-thin (3, 1) hypersurface, the (4, l)-dimensional DOD is negligible because of this thinness (see Fig. 
3). Other reasons are discussed below. 

Now we note that the Cauchy-Kovalevskaya theorem has no control whatsoever over these last two proper- 
ties for the analytic functions. There are furthermore two reasons why the analytic functions are undesirable. 
First, the analytic functions fail to be appropriate for any causal study since they are rigid (a change in a small 
region of an analytic data set induces an unwanted change of the entire data set). Without them one absolutely 
cannot have anything like the extremally general (and in particular signature-insensitive) Cauchy-Kovalevskaya 
theorem. Second, realistic^ and (or) interesting situations seldom involve analytic or smooth functions. For ex- 
ample astrophysical bodies have boundaries, thin matter sheets are discontinuous and the approach to spacetime 
singularities can be of low differentiability. 

Bruhat showed that the GR CP is well-posed in 1952 f|37j) and 1956 (|^), by making use of how the 
(n, 0; -1) EFE system (|15m7|) when cast in harmonic coordinates is of the correct quasilinear hyperbolic form 
£rA(2;^$^^VA$r)VrVA«'E = ^s(x,$A,VA$r) (where C^^ is a Lorentzian metric; both this and the function 
J-'s are smooth) to enable use of Leray's theorem |4H| I16j. which guarantees local existence and uniqueness, and 
furthermore continuous dependence on the initial data and the DOD property. 

The above proofs of the four well-posedness criteria for the GR CP are now usually done using Sobolev spaces 
jl6l 1141 1471 1151 1491 148j . These are appropriate because of the hyperbolic character of the Einstein evolution 
equations. At first they lend themselves to less involved proofs than Bruhat's, although if one seeks yet stronger 
results the functional analysis becomes extremally unpleasant. To recollect simply why Sobolev spaces are 
appropriate it suffices to consider a Klein-Gordon scalar in Minkowski spacetime. Given data on a bounded 
region 5*0 of a spacelike hypersurface Sq, one can draw the future DOD T>^{So) [Fig. 4a)] which is the region 
affected solely by this data due to the finite propagation speed of light. One can then consider the values of ip and 
its first derivatives on St = 'D~^{S()) r]T,t,t > 0. Then by Gauss's (divergence) theorem and energy-momentum 
conservation, 

/ TAct^t'' + / TAcl^'t^ = I TAct^f. (25) 

JS^^ JB JSt 

and the second term > provided that the matter obeys the dominant energy condition (DEC) 

— T^cu~' is a future-directed timelike or null vector V future-directed timelike u"^ (26) 

and that t"^ is timelike. Then the definition of the energy-momentum tensor gives 

/ ((^)2 + (VV')2 + mV) < / ((V')^ + (VV^)2 + mV) (27) 
J St J So 

Because each integrand is the sum of squares (which are necessarily positive), this means that control over the 
data on So gives control of the solution on St- The idea of a Sobolev norm is a generalization of these last two 
'energy' integrals [l^. For later use, the Sobolev class H'^ has a bounded norm of this type including up to wth. 
order derivatives. 

Whereas the GR CP comes with important well-posedness theorems jl61 1141 [TB] . the (r, 1; 1) problem is 
instead a complicated instance of a sideways Cauchy problem, a category of problem for which the possibility 
of such well-posedness theorems remains notably undeveloped.^*^ There is quite simply no established way to 

^Whereas Hawking and Ellis |14| argue that the choice of function space used to model nature does not matter since it is not 
experimentally-determinable and in any case is only an approximation due to QM, we would refine this argument to say that the 
difference between the analytic functions and other function spaces is important because of the rigidity problem. Beyond that, we do 
not know if the particular function spaces used to prove rigorous theorems about the EFE's and about the extendibility of spacetime 
|49| (of relevance to our Sec 5) may be substituted in these applications by approximations based on other function spaces. In this 
case, it may be a mere matter of convenience: one wants to use whichever sufficiently general function spaces accessibly give rise to 
the strongest possible theorems. 

^'^The few simple results known for fiat spacetime sideways wave equation problems are collected in 1501 : these results might serve 
as a starting point for the study of the much more complicated nonlinear sideways "GR CP" system. So little is known about 
ultrahyperbolic equations that we deem it not sensible to talk about (r, s; e) procedures at present for s > 1 or s = 1, e = — 1. In 
this paper we argue that even s = 1 e = 1 is hard enough! 



proceed. There is no known sideways (s = 1) analogue of Leray's theorem. There is no known notion of what 
function spaces might be appropriate for sideways Cauchy problems. We can explain however why the above 
introduction of Sobolev spaces does not carry over to sideways Cauchy problems. First there is no sideways 
notion of DOD to make the construction. Second even if we assume the higher-dimensional DEC holds, it would 
not give an inequality because the perpendicular vector in now spacelike [Fig. 4b)]. Third, we obtain a 
difference of squares rather than the sums in so the equivalent of the energy method's use of Sobolev norms 
is simply of no use to control the "evolution" given the data. 

3 Comparison of Methods of Data Construction 

Several methods have been proposed for the data construction step, including the elimination and conformal 
methods. The former methods are intuitive in that the prescribed quantities are all physical, but have a number 
of undesirable mathematical features. In contrast, Lichnerowicz's conformal method 36 is counterintuitive, 
since some of the prescribed quantities are unphysical, but exploits the mathematical properties of the constraint 
system in order to decouple it. Bruhat's argument ^38^ was in favour of the latter, however, some of her criticisms 
need to be elaborated upon to cover Magaard's elimination method j2ZI- Sec 3.1 covers this question and provides 
further criticisms in the s = 1 case. In Sec 3.2, we discuss the desirable features of York's development of the 
conformal method ^lE]) paying careful attention to the signature-dependent ones. The thin sandwich and 
Hamiltonian methods are briefly discussed in Sec 3.3. 

3.1 Elimination methods 

An example of this method of data construction is employed in the last part of the CM result ("Magaard's 
method" |271 I29j'). We discuss and compare this method with the standard method used for the GR IVP. 
Magaard's method (1963) treats the lower-dimensional metric hah as a known. Here we use a more transparent 
notation than Magaard's, in which the Gauss constraint ()15() takes the form 

G"^*^'^K,fei^,rf + ei? + 2p = (28) 

where G°^"^ = K'"}^'^ - W^h"'^ is the (undensitized) inverse of DeWitt's supermetric This is then split w.r.t 
some coordinate x\. Then K\\ can be isolated as the solution of the linear equation^^ 

IK^^G^^^^K^^ + G^^^^KruKx^ + G'''"^yK^^K^y = -(ei? + 2p) (29) 

(where u, x, y ^ 1), provided that it is possible to divide hy Z = G^^^^ Kuw Thus Ku is eliminated from the 
Codazzi constraint ()16() . which is then treated as a p.d.e system for unknowns U = {Ki^, some K^^ component 
denoted E}. On the face of it, this system satisfies the criteria for the CK theorem if one treats U' = {all the 
components of Kij bar Ku, Kiu and as known functions on xi and provided that the p.d.e's coefficients 
and the data are analytic. So a unique solution exists. The (generalized) CM result [231201 groups this and the 
local existence of a unique "evolution" to form the statement that a (r, s) spacetime with prescribed analytic 
metric hab may be embedded with an extra space or time dimension for any analytic functional form of the 
energy-momentum tensor. This statement suggests an incredibly rich collection of embeddings exist. However, 
we must point out that this is itself one of many difficulties associated with the CM result and its applicability, 
that we describe below. 

Bruhat had already considered the above method in 1956 ^] in the case of the GR IVP and pointed out a 
limitation on its validity. We find that her argument can be made for any signature or dimension: Z must be a 
function of at least one unknown, E, so we have no control over whether it is zero. Therefore we cannot guarantee 
the validity of the elimination procedure for Ku from the Gauss constraint. Although Magaard [2ZI finds a route 
round this problem, we find that this leads to two prices to pay later on . He starts with prescribed 'data for the 
data' i.e values of lA on some (n — l)-dimensional set Xi = {xi = 0}. The 'data for the data' is validly picked 
so that Z ^ Q on Xi, whereupon by continuity there is a thin region < < ry within which it is guaranteed 
that Z 7^ 0. Thus Magaard's procedure produces strips of data to be used in the embedding step. The first price 
to pay is that in general we cannot expect to be able to patch such strips together to make extended patches 
of data. For, since the strip construction ends where Z picks up a zero for some xi = r], while restarting the 
procedure with xi = r/ in place of xi = is valid, the two data strips thus produced will have a discontinuity 

^^The Gauss constraint is linear, not quadratic, in K\i by an obvious antisymmetry of the inverse supermetric. This spht uses 
normal coordinates so hia = and hn — 1. 



across xi = ij. So what one produces is a collection of strips belonging to different possible global data sets. 
The evolution of each of these strips would produce pieces of different higher-dimensional manifolds. So one is 
not in fact saying that an empty (n + l)-dimensional manifold surrounds any n-dimensional manifold. Rather 
one is saying that any n-dimensional manifold can be cut up in an infinite number of ways (choices of the xi 
coordinate) into many pieces, each of which can be separately bent in an infinite number of ways (corresponding 
to the freedom in choosing the components of Kab in W on each set of 'data for the data'), and for each of these 
bent regions we can locally find a surrounding (n + l)-dimensional manifold for every possible analytic function 
form of Tab (corresponding to the generalizability of the CM result). This apparent excess richness raises the 
question of physical significance of such a construction. 

Magaard's method does not state enough assumptions to make it rigorous. First, how far does 'data for the 
data' extend along xi = 0? Clearly the topology of the n-dimensional manifold is an important input, for if it 
is not compact without boundary, there is a missing boundary or asymptotics prescription required. Also the 
topology of the xi = set itself has not been brought into consideration (for example continuity is not guaranteed 
were this to contain loops). 

Secondly, we specifically consider the (n, 0; -1) and (r, 1; 1) problems as separate cases, since we find that 
there are implicit ways in which signature plays a crucial role even for the CM result. For (n, 0; -1) one might 
worry that the data construction for a strip is flawed because the data problem in question is elliptic and hence 
naively requires a global treatment. However, we are able to dispel this worry once we treat York's method 
below. For (r, 1; 1), the use of the strip < xi < r/ as "evolutionary" data is generally invalidated by the 
information leak construction in Fig. 5, unless one has had the luck to construct a full global data set. This 
would, however, require the building of the data encountering 1) no zeros of Z (which is thus the second price 
to pay), 2) no asymptotic problems. Thus in general there is a severe problem with building (4, 1) spacetimes 
from (3, 1) ones using the CM method. 

Thirdly, Magaard's method lacks any n-dimensional general covariance since it involves the choice of a 
coordinate xi and the elimination of a 11-component of a tensor. 



3.2 Conformal methods 

The arguments of Lichnerowicz and Bruhat [HH] led to the GR IVP taking a very different route from 
the above sort of brute-force elimination methods. The conformal method of Lichnerowicz and York is instead 
preferred. This leads us to ask to what extent conformal methods can be adapted both to heuristic and to general 
constructions for e.g (r, 1; e) data. Therefore our treatment below is as far as possible for the general (r, s; e) 
case. 

In the conformal method, one chooses to treat hij as a known metric which is moreover not the physical 
metric but rather only conformally-related to it by 

kj = (j^'^h^j. (30) 

We work in terms of KJ^ and K, and permit K'^^^ , f and p to conformally transform according to 

i^Tu- = /-^'^KTii ^ f = p=^-p, (31) 

whilst crucially demanding the constant mean curvature (CMC) condition 

K = hypersurface constant (32) 

holds and is conformally-invariant. One then demands that the (raised) Codazzi constraint H20|) is to be 
conformally-invariant. Since 



one requires that 



(33) 



n -I- ? 

=C-2r? = 6 (34) 



Furthermore one demands that the conformally-transformed Gauss equation 

R-r]{n-l)^ +7]{n-l)(l- r?^^) + e(M02C-^ - p^^'n + 2p,/)'^+'') = (35) 



(where M = K^oK'^ and /i is some "York extrinsic dynamical variable" proportional to K) is simplified by being 
made to contain no \D(f>\'^ term, so that rj = ( = —2 and ^ = —2^^ [by making use of (|3i|) ]. Now, regardless 
of (r, s; e), provided that the (n + l)-dimensional DEC is to be preserved by the conformal transformation, p'^ 
must conformally-transform like j"ja, implying oj = —2^^. Then ()35() becomes the (r, s; e) version^^ of the 
Lichnerowicz equation 

o e n — 2 / A n-\ r, 4 n n-l \ 

= — ^(p [-eR - M(j)~^~ + fi^(t>~ - 2p(f>~^—j (36) 

for the conformal factor (p. We observe that a number of the very attractive features of the York method are 
(r, s; e)-independent. First, it decouples the solution of the Codazzi constraint from the solution of the Gauss 
constraint. The former proceeds by a traceless-transverse (TT)-traceless-longitudinal (TL) splitting jl8j 

(for some vector potential W^), which along with the trace-tracefree split is conformally-invariant and thus 
unaffected by the solution of the latter, which has become the p.d.e (|36)) for the conformal factor. The simplest 
case IS 

j^Tij ^ K™^, f = which at most requires solution of a first-order linear equation. 

Second, the choice in using the scale-scalefree decomposition of the metric and the trace-tracefree and TT- 
TL decompositions of Kij are all decompositions into irreducibles and thus n-dimensionally generally-covariant 
choices, a decided advantage over the Magaard method. 

Third, one can attempt to preserve the CMC condition away from the 'initial' hypersurface Tq, which can 
be done provided that the equation 

D^a + a ( eK"^ — R — ^"'^ + S \ _ ]^ypgj,g^j,£g^gg constant (38) 
V n-l J 

[from the first form of H18|) ] can be solved for the lapse a. This choice of lapse is the CMC slicing gauge. This 
deliberately defocusses geodesies, thus by definition preventing the breakdown of the coordinate system due to 
caustic formation |361 143j . and so enhancing the practical longevity of the evolution. However, York's method 
is not absolutely general, for some spacetimes have no CMC slice to identify with Tq in the first place, and in 
others the CMC slicing cannot be maintained to cover the whole spacetime. Furthermore, these results depend 
on the asymptotics assumed. ^'^ However, in the usual (3, 0; -1) case, this method (and its variants) is widely 
accepted as a practical method by the numerical relativity community, for example in the study of colliding 
compact astrophysical objects |52( 1201 153j . Of importance for the viability of this method in this application, 
recent evidence suggests that problems involving gravity wave emission can be treated thus j54j . 

However, some of the other attractive features used to gain control in the York method are only known to occur 
for the s = case and have mainly been studied in depth for the compact without boundary and asymptotically- 
flat cases of (3, 0) data construction. First, the study of the usual (3, 0; -1) Lichnerowicz equation has been based 
on its ellipticity (and elliptic methods are absolutely not generalizable to the 'hyperbolic' equations^^). The (3, 
0; -1) Lichnerowicz equation is well-studied, including with most fundamental matter flelds 15^, and for Sobolev 
spaces matching those then used in the GR CP 56, 15^. Second, the study of the CMC slicing equation has been 
based on its ellipticity, so the theory of existence of CMC slicings [ISl EH is likely to be signature-dependent. 
Third, the method depends on the asymptotics used and has only been studied in the (3, 0) compact without 
boundary |581ll5j and asymptotically fiat j591ll5j cases. Fourth, for s = there is the useful property that certain 
local data patches can be proved to suffice for the treatment of an astrophysical problem, by building on the 
notion of DOD. As mentioned, this sort of technique is also applicable to protect pieces of local data obtained 
by Magaard's method on a (r, 0) hypersurface. 

Investigation of whether the s = 1 'wave Lichnerowicz equation' has good existence and uniqueness properties 
could be interesting. The natural setting for this is as a Cauchy problem (for even the 2-dimensional wave equation 
is ill-posed as a Dirichlet problem). Assuming that there exists a Cauchy surface in the (3, 1) spacetime sense, 
one can attempt forward and backward evolution to produce a global data set (assuming also that the decoupled 
procedure for finding /Cp^ also yields a global solution). One must remember however that the next stage is still 
to be a sideways Cauchy problem. Whilst this suggested procedure for the data could conceivably produce global 

^^The main feature due to s is hidden in . In the specific case of s = 0, is the elliptic Laplacian operator denoted by A. 
^^It was these caveats, along with the reliance on s = of York's method and restrictions on the values R can take for some 
subcases of the Lichnerowicz equation, that made us suspicious of the supposed generality of the Magaard method. 
^^Unlike ellipticity, hyperbolicity is actually difficult to pin down for nonlinear systems |45|. 



data sets in some subcases and thus avoid the information leak problem, the other difficulties we described in 
Sec 2.4 remain. Therefore in our view the 'wave Lichnerowicz equation' is unlikely to be suitable as a general 
method. We favour instead the (4, 0; -1) approach in Sec 6. That said, the 'wave Lichnerowicz equation' may 
still serve as the basis of a useful heuristic method. We consider this in Sec 4.7. 



3.3 Other methods and formulations 

We finally discuss the thin sandwich method and why we do not provide an (r, s; e) Hamiltonian formulation. 

In GR the lapse may be algebraically eliminated from the action by use of its own variational equation to form 
a reduced Baierlein-Sharp- Wheeler action [301 • Suppose that {hij, hij) is prescribed on a hypersurface along 
with p and ji. The thin sandwich conjecture [39j is that in reasonably general circumstances there exists a unique 
solution to solving the variational equation for the shift (this equation replaces the Codazzi constraint), as a 
differential equation for the shift itself. Once this is done, the lapse may be reconstructed from the equation used 
to eliminate it, and the extrinsic curvature then follows from its definition. The lapse may then be reconstructed 
from the equation used to eliminate it, and the extrinsic curvature then follows from its definition 

We note that this procedure is (r, s; e)-independent. Thus we pose the conjecture for general (r, s; e). This is 
also considered for thin matter sheets in Sec 4.7. To date we are only aware of (3, 0; -1) papers j31j treating this 
conjecture, which contain some favourable and unfavourable results obtained mostly by use of elliptic methods. 
In fact, there are two variants of the conjecture: i) the thick sandwich conjecture in which lower-dimensional 
metrics are prescribed on two nearby hyper surf aces, ii) The thin sandwich conjecture proper, as stated above. 

A (3, 1; 1) Hamiltonian formulation is not provided because although it is insensitive to e j32|, it is sensitive 
to s through the intimate involvement of the space of geometries. Whilst the usual Hamiltonian formulation is 
based on the space of Riemannian geometries (superspace) j39( I51j. the space of semi-Riemannian geometries is 
reported to be not even Hausdorff (60) . 



4 Application of (r , 1 ; 1) Methods to Thin Matter Sheets 

First, we present a simple instance common in the higher-dimensional literature to which Sec 2.1-3 is applicable. 
The warpfactor split (7| 

,co = ( f »-(.n,,);,,(,n) ) . W(.-,0) = 1 (39) 

is a simple subcase of the z-dynamics scheme, in which the metric is allowed to z-evolve only in its scale, away 
from the z = hypersurface where it is taken to be known. Then leads to df^^'^/dz = — /^/^<1>/C which gives 
the equation 

^ = -*/C (40) 

oz 

for the warpfactor. For example, using the ansatz /CpA = CfrA in (|4U|) gives an exponential Randall-Sundrum 
type warpfactor. Whereas the split (|39|) does not cover very many cases, our scheme exhibits generalizations for 
it: to permit the whole metric to evolve and to recognize the gauge freedom in /3r, which should ideally be used 
to separate coordinate effects from true physics in the spirit of |i3] . 

Full, overtly z-dynamical schemes are used in particular examples for domain walls j34j . braneworld black 
holes (such as for the pancake or cigar bulk horizon shape problem j21j ) and for braneworld stars |221 I23j . 
Whereas in the Randall-Sundrum model the higher and lower-dimensional cosmological constants balance out 
leaving vacuum (Minkowski) on the brane, more generally a brane would consist of a thin sheet of matter - a 
junction. We study such [Shiromizu-Maeda-Sasaki (SMS)-type] braneworlds below, ""^^ starting first however with 
a careful recollection of where the underlying Israel junction conditions come from. 



4.1 Origin of the Junction Conditions 

Assume that we have a {D - l)-dimensional thin matter sheet in a D-dimensional bulk. In all cases considered, 
the bulk's extra dimension is spatial (e = 1). Our discussion follows ^S] most closely whilst keeping the unraised 
index positions of Israel's original work jl2j . 

^^This section presents the material announced in |70| in a more advanced form and provides the derivations of the results along 
with a number of additional points and further new results. 



Whereas the requirement of well-defined geometry dictates that the metric is continuous across the thin 
matter sheet yielding the junction condition (j.c) 



[/rA]l = /r+A - /Fa = 0, (41) 

discontinuities in certain derivatives of the metric are permissible. Consider then the 3 projections of the Einstein 
tensor Gab- We use the e = 1 cases of the Codazzi and Gauss constraints (fT^ and ((T^ for Gj-± and G±± 
respectively. For GrA, the following construction is used. 

One begins by writing down the contracted Gauss equation ((TJ and subtracts off ^/rA times the doubly- 
contracted Gauss equation - ((TU)) : 



T^rA - T^±r±A = Rta - K.1Cta + ^aa 

2 



-f/rA + H^Uta -l/rA + ^^fvA. 



(42) 



The following steps are then applied. 

Step 1: The Ricci equation © is used to remove all the 7?._Lr±A- 

Step 2: The contracted Ricci equation (jUJ is used to remove all the 1Z±±. Thus 



7^^A 

2 



l/rA + 



h /Cr /Cah 



_ ( Rta \ - /C/CrA + ^r^A^AA ..„n 
fvA V "2/rA ; + — 2 — /rA. 



This is then rearranged to form the "GR CP" geometrical identity 

^ K? + jColC Sp^TA - DrD/^a - {S^K. - D'^a)frA 

^rA = GrA - ^^rA + 2/Cr /Caa H /rA H • (44) 

2 a 



lim 

Then performing / GAsdz one obtains the j.c's 

C ^ O J — e 



[G±±]t = 0, [Gr±]t = 0, (45) 
[GrA]t = [ICta - fvAK]t. (46) 

The derivation of this last equation makes use of normal coordinates (in which case the hypersurface derivative 

5^ becomes the normal derivative 

p az I 



S3 becomes the normal derivative ^) and the rearrangement to the 'Israel' geometrical identity 



GrA = GrA + /C/CrA + 2/Cr'^/CAA H /rA + tt (^fa - JtaK) (47) 

2 oz 

via the definition of extrinsic curvature (|2j to form the complete normal derivative ^(/CrA — fvA^) ■ 

Step 3: One then further assumes that the (4, l)-dimensional EFE's GvA = Tta hold. If one then additionally 

lim 

uses the thm matter sheet energy-momentum yAB= T^sdz. one obtains the j.c's 

= , = 3^r± , (48) 

y 

L> - 2' 

(performing a trace-reversal to obtain the last equation). 



[ICvAt =[yTA- 7^/rA ) (49) 



4.2 SMS's Braneworld EFE's 

We next recall the method SMS use to obtain their BEFE's 0. They begin by forming the (3, l)-dimensional 
Einstein tensor GrA just like we obtain ()42() above. SMS then apply three steps to this equation. 
Step SI: Using the definition of the Weyl tensor, 7^_Lr±A is replaced by the electric part of the Weyl tensor, 
£ta = C_Lr±A and extra terms built from the projections of IZab- 

Step S2 ( = Step 3 of the above subsection): The (4, l)-dimensional EFE's are then assumed, which permits 
one to exchange all remaining projections of TZab for (4, l)-dimensional energy-momentum terms. Only when 
this is carried out does ()42|) become a system of field equations rather than of geometrical identities. We refer to 



the field equations at this stage as "timehke hypersurface EFE's" (THEFE's)/^ as opposed to the braneworld 
EFE's which arise at the next stage. 

Step S3: A special subcase of THEFE's are braneworld EFE's (BEFE's), which are obtained in normal coordinates 
by choosing the (thin) braneworld energy-momentum tensor ansatz 

Tab = yAB^z) - Kqab , yAB = {Tab - X/ab) , Tabz^ = 0, (50) 

where T^a is the energy-momentum of the matter confined to the brane. This is a specialization due to the 
specific presence of (3, 1) and (4, 1) cosmological constants A and A, and by A being the only bulk contribution. 
A more precise formulation of this, and generalizations, are the subject of Sec 4.5. 

One then adopts the j.c's (|^T|) . and with the additional supposition of Z2 symmetry^'': 



- ^rA = ^rA = "''^rA =^ ^rA = --^ (yrA - y/rA^ = --^ (tta 5 — /rA ) , (51) 

where the 5-dimensional gravitational constant k| has been made explicit. Then SMS's BEFE's read 

GrA = -^^rA^ + QrA^ ~ "^rA, (52) 

where QrA{T) and Lya{T) are the terms quadratic in, and linear together with zeroth order in T respectively, 
given by 



^SMS _ 4 

VrA — 



^TrA - ^TrnT^A + ( ^ - ^ ) fvA 



(53) 



Lff = -f (^A + j /rA + f ATrA. (54) 

As opposed to the (3, l)-dimensional EFE's, SMS's BEFE's are not closed since they contain the unspecified 
electric part of the Weyl tensor SrA- Although it also contains 15 equations, the SMS BEFE-Gauss-Codazzi 
system^^ is not equivalent to the (4, l)-dimensional EFE's: indeed SMS write down further third-order equations 
for the "evolution" away from the timelike brane of Sta, by use of the z-derivative of the contracted Gauss 
equation ((Tj), Bianchi identities and the Ricci equation @. This then involves the magnetic part^^ of the Weyl 
tensor ;BAAr = C^rAA, the "evolution" of which follows from further Bianchi identities. This full brane-bulk 
SMS system is then closed. 

Step 4: In practice, however, instead of the difficult treatment of this third-order system, other practitioners 
have often worked on the SMS BEFE's alone. This involves either the ad hoc prescription of the functional form 
of SrA (sometimes taken to be zero^'^ ). In fact it is often first decomposed according to a standard procedure 
p'2l I63j . Because the original functional form is unknown, the functional forms of each of the parts defined by 
the decomposition is also unknown. Some of these parts are set equal to zero whereas other parts are taken to 
have other functional forms (in particular a radiation fluid term). These terms are then argued to be small in 
the circumstances arising in the inflationary [H^j and perturbative [SSI treatments. However other interpretations 
may be possible. Having dealt with S^a in one of the above ways, the form (53) of Qp^^ then often taken to 
be uniquely deflned and the starting-point of many works on brane cosmology |M1 1^ 1^ 1^ . 

However, SMS's procedure is far from unique. It turns out that there are many reformulations of the BEFE 
arising from geometrical identities. Each has a distinct split of the non-Einsteinian BEFE terms into 'bulk' and 
'brane' terms. Whereas all these formulations are clearly equivalent, their use helps clarify how to interpret 
SMS's braneworld. Were one to truncate the 'bulk' terms in each case (in direct analogy with the usual practice 
of throwing away the Weyl term), then the BEFE's obtained in each case would generally be inequivalent. 

^^This choice of name reflects their superficial resemblance to the (3, f )-dimensionaI EFE's, although as discussed below, this 
resemblance is often only superficial, including in the SMS case. 

^^The difference in the sign of 1491 between this paper and SMS's paper is due to our use of the opposite sign convention in the 
definition of extrinsic curvature. We compensate for this in subsequent formulae by also defining Kab = ~^^t, rather than +K^^. 

^*The 5 Gauss-Codazzi equations may be seen as consistency conditions on the brane. For jr = the 4 Codazzi equations 
correspond to energy-momentum conservation on the brane. 

^^As the 5-dimensional alternating tensor has 5 indices, one has two choices for the number of indices in what is to be taken to be 
the definition of the magnetic Weyl tensor. The above is the 3-indexed definition; the other possible deflnition is //^"^ — e^'^^^^Brsn- 

^"it is zero for example when one makes the restrictive presupposition of a conformally-flat bulk, such as in AdS spacetime. 

■^^For example the 'Weyl charge' for black holes in 



4.3 Ambiguity in the Formulation of the BEFE's 



The Weyl term in SMS's BEFE's has been the subject of much mystery. How should it be interpreted? Is it 
right to throw it away and if so under which circumstances? We emphasize that our stance is broader than 
merely about what functional form is allocated to the Weyl term (e.g whether it is zero everywhere). It is about 
how formulations can be chosen in which the BEFE's are not explicitly formulated in terms of a Weyl term. 
We first remove some misconceptions as to how SMS's procedure leads to a BEFE containing a Weyl term. 
Does it have anything to do with the modelling of braneworld scenarios? No, for the Weyl term is already in 
the SMS THEFE before the braneworld energy-momentum ansatz is invoked. Furthermore, all the procedures 
used in SMS's method are independent of signature and dimension. Thus this issue of a Weyl term must have 
already arisen long ago in the study of the GR CP. So why is there no manifest Weyl term piece in the GR CP 
formulation of the EFE's? The answer is simple. Indeed, we have already seen the answer in Sec 4.1 since it is 
fundamentally tied to how the crucial junction condition l|49j) is obtained in the first place!. In the "GR CP" 
and Israel procedures, one uses the Ricci equation ^ to remove the 7^_Lr±A term. If there is no early use of the 
Ricci equation, one is left with S-pA in the THEFE's, which requires later use of the Ricci equation to "evolve" 
it. 

So there is a choice as to whether one formulates the BEFE's with or without an explicit Weyl term. In the 
usual treatment of the split of the EFE's (Sec 2.2), one does not use an explicit Weyl term. In the (4, l)-d case, 
this gives a well-understood system of 15 p.d.e's in the variables /fa- The option of using an explicit Weyl term 
gives a considerably larger, more complicated system of p.d.e's with variables frA, ^VA and iSpAA- 

In fact a similar scheme known as the threading formulation 68_ is sometimes used in the usual (3, 0; -1) 
application. The idea behind this formulation is to treat as primary the geodesic congruences perpendicular to 
the foliation rather than the foliation itself. One then uses only that information on the hypersurfaces that arrives 
along the impinging geodesic congruence, on the grounds that this is the physically-significant information. Thus 
it is a 'deliberately incomplete system' from the foliation perspective. Whereas this and the SMS formulation are 
similar in their use of Weyl variables, the SMS system does not appear to be a threading formulation. Also, in 
any case the idea of having a deliberately partial system in the threading formulation is clearly tied to signature- 
specific physical reasons which do not carry over to the signature relevant to SMS's equations. We also note that 
some other third-order reformulations of the (3, 0; —1) split of the EFE's are sometimes used to seek to cast the 
EFE's into hyperbolic forms that manifestly have theorems associated with them |45| . Whereas this is precisely 
the sort of result that is spoiled by considering instead a sideways split^'^, it serves to illustrate that what at 
first seems a 'mere reformulation' of a set of equations can in fact be used to prove highly nontrivial theorems. 
So, whereas similar complicated formulations have been used elsewhere in the GR literature, SMS's unstated 
motivation to have a complicated formulation does not coincide with the motivation elsewhere in the literature. 
Below we bring attention to many reformulations of SMS's system, so we ask: what is the motivation for the 
original SMS formulation? Should the use of some simpler second-order formulation be preferred? Is SMS's 
formulation or any other third- or second-order formulation singled out by good behaviour, either in general or 
for some particular application? 

Also, from first principles the SMS procedure to obtain their BEFE is quite complicated. For, since they 
use the j.c obtained by the Israel procedure, their procedure actually entails beginning with the whole Israel 
procedure (Steps 1 to 3 of Sec 4.1), and then choosing to reintroduce TZr±A± and TZ±± by reverse application 
of Steps 1 and 2. This is followed by the Weyl rearrangement (Step SI), the use of the EFE's (Step S2) and the 
substitution of the j.c into the extrinsic curvature terms in the braneworld ansatz (Step S3). However, despite 
being complicated, all is well with SMS's scheme since any BEFE's obtained by other such combinations of 
careful procedures will always be equivalent because the different steps are related by geometrical identities. 

Step 4 however is not an instance of being careful as it is a truncation. Our first point is that whereas in 
SMS's formulation the non-Einsteinian terms in the BEFE might be regarded as a bulk-like £rA and a term 
quadratic in the brane energy-momentum, in other formulations the content of these two terms can be mixed up. 
In general, BEFE's contain a group of non-Einsteinian 'bulk' terms we denote by B-pA (which include both Weyl 
terms and normal derivatives of objects such as the extrinsic curvature), and a group of non-Einsteinian 'brane' 
terms that depend on the brane energy-momentum. Thus any temptation to discard the Weyl term in the SMS 
formulation (on the grounds that it involves the unknown bulk over which one has no control) should be seen in 
the light that if one considered instead a reformulation, then there would be a similar temptation to discard the 
corresponding 'bulk' term, which would generally lead to something other than the Weyl term being discarded. 

■^^Having argued that the second-order sideways "GR CP" is not known to be well-posed, we should add that it is unlikely that 
third-order formulations are likely to be better-behaved in this respect. 



Thus for each formulation, the corresponding truncation of the 'bulk term' would result in inequivalent residual 
'braneworld physics'. This is because there are geometrical identities that relate 'bulk' and 'brane' terms, so 
that the splits mentioned above are highly non-unique and thus not true splits at all. We take this as a clear 
indication that any such truncations should be avoided in general. Instead, the full system must be studied. 

Our second point is that each possible bulk spacetime may contain some hypersurface on which a given B^a 
vanishes. Then if one identifies this hypersurface with the position of the brane, one has a solution of the full 
brane-bulk system and not a truncation. For example, in any conformally-flat bulk, by definition Cabcd = 
and therefore i?p^^ = —SrA = on all hyper surf aces. Thus any of these could be identified as a brane to 
form a genuine (rather than truncated) — <?rA = braneworld. From this, we can see that the SMS formulation 
is particularly well adapted for the study of conformally-flat bulks such as pure AdS. This motivates SMS's 
formulation as regards this common application. However, also consider repeating the above procedure with 
some -BrA 7^ —^VA- This would correspond to a genuine (rather than truncated) braneworld model with distinct 
braneworld physics from that given by SMS's particular quadratic term. Note that given a model with some 
-Br A = 0, the BEFE formulation for which -Br A is the bulk term is particularly well adapted for the study of that 
model. Thus different formulations may facilitate the study of braneworlds with different braneworld physics. 
In the context of conformally-flat spacetimes, it is probably true that the —Sta = braneworlds outnumber the 
braneworlds for which any other (or even all other) -BrA = since these other conditions are not automatically 
satisfied on all embedded hypersurfaces. Rather, each of these other conditions constitutes a difficult geometrical 
problem, somewhat reminiscent of the question of which spacetimes contain a maximal (K = 0) or CMC slice 
j69j . However, generic spacetimes are not conformally-flat. For a generic spacetime, we see no difference between 
the status of the condition '— <SrA = on some hypersurface' and the condition 'any other particular -BrA = on 
some hypersurface'. Because braneworlds constructed in each of these cases have a different residual quadratic 
term and thus a propensity to have distinct braneworld physics, and because we do not know how frequently each 
of these cases occur, we question whether anything inferred from conformally-flat models with the SMS quadratic 
term need be typical of the full SMS brane-bulk system. Confirmation of this would require study of the full 
range of difficult geometrical problems -Br a = 0, and the construction of concrete examples of non-SMS quadratic 
term braneworld models together with the assessment of whether their braneworld physics is conceptually and 
observationally acceptable. 

For the moment we study what is the available range of reformulations and thus of -Br A- To convince the 
reader that such reformulations exist, we provide a first example before listing all the steps which are available 
for reformulating the BEFE's. 

Assume we do not perform all the steps implicit within SMS's work but rather just the Israel steps to obtain 
the j.c and then use it in the field equation (|T7|) that gave rise to it (as done in ^01)) or (as done below) use it 
in the "GR CP" field equation following from H44() . In other words, why not apply the braneworld ansatz to e.g 
the Israel or "GR CP" formulations rather than to the SMS formulation? In the "GR CP" case we then obtain 

GrA = LrA + Qfa + -BrA , with (55) 



QrA = -^ [seTr^TnA - ISTTta + (9T o T + r2)/rAj (56) 
LvA = ^iT\- 2\^)fTA + ^^[TvA - (A + A)/rA] (57) 

BrA = /r.^-^. (58) 

oz oz 

This example serves to illustrate that choosing to use a different formulation can cause the 'brane' quadratic 
term QrA{T) to be different. Also note that this formulation makes no explicit use of the Weyl term. Thus this 
BEFE, along with the Gauss and Codazzi constraints, forms a small second-order system, in contrast with the 
much larger third-order SMS system. 

Now we further study the list of steps [70] which may be applied in the construction of BEFE's. 
Steps SI and 3 together mean that the Weyl 'bulk' term E^a is equivalent to the Riemann 'bulk' term together 
with matter terms. This swap by itself involves no terms which are quadratic in the extrinsic curvature. 
Step 1 says that the Riemann 'bulk' term is equivalent to the hypersurface derivative of the extrinsic curvature 
together with a /Crn^^A term. 

Steps SI and 3 together say that the hypersurface derivative of the trace of the extrinsic curvature is equivalent 
to a matter term together with a /C o /C term. 

Furthermore, one can use both Steps 2-3 and Step 1, on arbitrary proportions (parametrized by ji and v) of 
'^±r±A and of 7^_l_l: 



GrA = QvA - (1 + i^)^±r±A+(l - /u)7^x±/^A + ^ ^{S^ICva - DrD^a) + fi{6^IC - D^a)frA 



+/C/CrA + {1^- l)/Cr"/CAn - ^fvA + (h - o /C/fa 



(59) 



This introduces freedom in the coefficients of the /Crn^'^A and /Co/C contributions to the quadratic term QrA(^) 
in the THEFE's. We next find further freedom in Qta by choice of the objects to be regarded as primary. 

We are free to choose a 'bulk' term described by hypersurface derivatives 6^ (which are partial derivatives ^ 
in normal coordinates) of objects related to the extrinsic curvature /CrA by use of the metric tensor (including 
its inverse and determinant /). The underlying reason for doing this is that it is just as natural to treat such an 
object, rather than the extrinsic curvature itself, as primary (see below for examples). 

Upon careful consideration, there are three separate ways such objects can be related to the extrinsic cur- 
vature: raising indices, removing a portion of the trace by defining ,,/CrA = ^rA — V^frA, and densitizing by 
defining ^/CrA = ^ ^VA- The hypersurface derivatives of these objects are related to those of the extrinsic 
curvature by 

SplCvA = S^ifruK^A) = fru^^lC^A - 2a/Crn^"A, (60) 
S^ICrA = ^p^TA + ^ {^^^frA - 2a/C/CrA) , (61) 

S0lCrA = Ur^'S^if^'ICrA) + 2aiilClCvA. (62) 
Further useful equations arise from the traces of these: 

5^K: = f^S^KvA + 2a/C o /C, (63) 

= if^^0vf^rA - 2aviJC' - /C o /C)] , ,? / ^ (64) 

S^fC = {fr^H^U^^K,) + 2ai2lC\ (65) 

where the in (|BH) and has been obtained via 

The following examples of |/CrA illustrate that the use of such objects is entirely natural: \iJCya is the /Cp^ 

1/2 

commonly used in the IVP literature, and g /CrA appears in the guise of forming the complete normal derivative 

1/2 

in the Israel procedure. Also, the "gravitational momenta" are pya = — i ^VA- 

The above thorough consideration of possible 'bulk' terms permits all four THEFE terms homogeneously 
quadratic in the extrinsic curvature to be changed independently. One may think that we have a redundancy 
in providing 8 ways to change only 4 coefficients. However, one can afford then to lose some of the freedoms by 
making extra demands, of which we now provide four examples of relevance to this paper. First, one could further 
demand that there is no Weyl term in the THEFE's (as discussed in Sec. 4.4). Second, unequal densitization 
of /CrA and /C (^i 7^ ^2) corresponds to interpreting the fundamental variable to be some densitization of the 
metric rather than the metric itself. Whereas this is again a common practice (for example the scale-free metric 
of the IVP literature is /~"/rA in dimension n), the use of such an object as fundamental variable does appear 
to complicate the isolation of the (3, 1) Einstein tensor truly corresponding to this fundamental variable. Thus 
this option is not pursued in this paper. Third, one may start by declaring that one is to use particular well- 
known primary objects (such as the "gravitational momenta") and still desire to be left with much freedom of 
formulation. Fourth, one could declare that one is to use the raised objects given by H63|). in which case the 
further ability to change coefficients by use of (|61j) is lost, since moving a Kronecker delta rather than a metric 
through the derivative clearly generates no terms quadratic in the extrinsic curvature. 



4.4 Examples of Formulations of BEFE's with no Quadratic Terms 

As a consequence of the above freedoms, there are many formulations in which all four coefficients vanish, and 
hence (5rA(/C) = 0. From this it follows that QrA{T) is zero [and it is easy to show that all instances of 
Qta{T) = follow from QrA(/C) = 0]. Thus it suffices to seek for cases of THEFE's with QrA(/C) = to obtain 
all cases of BEFE formulations that have no quadratic term. We now motivate these formulations and then 
choose to exhibit three that comply with some of the extra demands in the previous paragraph. 

The diversity of 'brane'-'bulk' splits ensures that truncations such as Step 4 produce all possible combinations 
of quadratic terms as residual 'braneworld physics'. Alternatively, one may suspect that there might be solutions 
to the full brane-bulk system that just happen to have a particular -BrA = on some hypersurface which is then 
identified as a brane. We speculate that each of these situations will often lead to different answers to questions 
of physical interest. Whereas most Friedmann-Lemaitre-Robertson- Walker (FLRW) perfect fluid models with 



equation of state P = (7 — 1)/) arising thus will be similar, differences will be more salient in models with more 
complicated equations of state, in perturbations about FLRW (as started in [5311^1721 in the SMS-adapted case) 
and in anisotropic models (as started by [HZl [Z3 in the SMS-adapted case). These in turn constitute natural 
frameworks to seriously justify the late-time emergence of FLRW behaviour and the likelihood of inflation |6t)j 
as well as the study of singularities 73, 72 on the brane. We emphasize that, for a satisfactory study of whether 
any particular full brane-bulk (as opposed to truncated) case leads to any differences from the hitherto-studied 
—SrA = case, one would require a full brane-bulk solution explicitly constructed to satisfy some i?rA = 
on some hypersurface within a particular given bulk spacetime. Since we currently have no such example, our 
arguments currently only support the far simpler idea that truncation should be avoided. 

We illustrate that in different formulations, the truncation of the corresponding 'bulk' terms can lead to 
big differences in the residual 'braneworld physics', without any of the above lengthy calculations. We do this 
by formulating the 'bulk' part so that there is no corresponding Q term at all. Thus these truncations give 
the 'p' of standard FLRW cosmology rather than the '/) and /j^' of brane cosmology |64 | I62[ l63] . As a result 
whether we have a '/? and /3^' brane cosmology depends on the choice of formulation. So we argue that since the 
SMS procedure followed by truncating the Weyl term is a hitherto unaccounted-for choice out of many possible 
procedures, then adopting the particular homogeneous quadratic term of SMS (often taken as the starting-point 
of brane cosmology) appears to be unjustified. Rather, we conclude that no particular truncation should be 
privileged as the act of truncation imprints undesirable arbitrariness into the study of the truncated system. 
Whereas the (3, l)-dimensional trace of if^rA happens to be zero (including use of antisymmetry) and thus 
might^^ phenomenologically look like pure radiation fluid to observers on the brane, other bulk characterizations 
would typically not look like a pure radiative fluid. This may open up phenomenological possibilities. 

Also, before further study of SMS's full (untruncated, third-order) system is undertaken, some of the re- 
formulations along the lines suggested in Sec 4.4 might turn out to be more tractable. In particular those 
reformulations which fully eliminate 7^_Lr±A by the early use of the Ricci equation are already closed as second- 
order systems. These include the Israel formulation in (70], the "GR CP" formulation (|55H58|) . and our second 
and third examples below, which contain neither a Weyl term nor a quadratic term. 

For our first example, we take as the primary object the antidensitized extrinsic curvature K^^i, = so that 
the 'bulk' term is (partly) a combination of this object's normal derivatives. The corresponding BEFE's are: 

GrA = Cta + Br A (66) 

^rA = ^ + ^ (5T^x - T) /rA , Br a = -2£rA + // - ^/"^^/pa) (67) 

(where we have chosen to remove all projections of TZab by the EFE's). 

To derive this, take (|59|) in normal coordinates. Choose to convert all of the ^ into /^^ ^^z^ by l|63|) : 

(68) 



GrA = ^rA - (1 + ^^)7^±r±A + (1 - Ai)^±±/rA + [i^^ + ^H/rA 
+/C/CrA + {1^- l)/Crn/CAn - ^fvA + + m)/C o /C/rA 

Now choosing the primary object to be some densitized ^i^rA by H62jl we have 

GrA = ^rA - (1 + i^)n±T±A + (1 - M)^±±/rA + ^ {^^^ + ^/^^^%^^ 

+(1 + 2u01ClCrA + (i^ - l)/Crn/CAn + (2//^ - ^)/CVrA + {\ + f^)JC o tfrA, 



(69) 



so clearly fi = u = 1 and the antidensitization choice of weighting C = — ^ ensure that QrA(^) = 

The following examples arose from asking if it is possible to find examples in which neither QrA{T) nor £rA 
feature. We found the following BEFE's: 

GrA = ^rA + Br A (70) 
LrA = 'rrA + l (tx± - I) /rA , ^rA = ^ (^Ij^^a - ^/ha) (71) 

by considering as our primary object the densitized extrinsic curvature with one index raised, fC°'b = VhlC"'b- 
We obtained these BEFE's by arguing as follows. In order for the BEFE's to contain no Weyl term, is fixed 

■^■^This is if one treats the geometric content of the Weyl tensor as an effective or 'induced' energy-momentum, and furthermore 
assumes that this ought to behave hke a perfect fluid. Maartens' decomposition of £rA |63) . based on the first assumption, permits 
the removal of the second assumption. 



to be —1. Then the only control over /Crn^'^A is from raising by (|6U|) . It is easy to show that this raising must 
be applied to the whole ^^j^ in order for the coefficient of /Crn^^^A to be zero. Then using jj/C^a does not 
change any terms quadratic in the extrinsic curvature. Also, use of distinct densities for /C and /C^a does not 
appear to make sense since both quantities are related to /CrA by a single use of the inverse metric. Although all 
these restrictions make the outcome unlikely, the use of (|62|) and (|H5|) alone suffices to obtain the above example: 



GrA 



= Gta + (1 - /^)7^±±/^A + (-%^ + /"f /rA 
+(1 - 20/C/CrA + (2e/x - ^)/C2/rA + - /u)/C o /C/rA 



(72) 



which has no quadratic terms if fi = 1/2 and = 1/2 ('densitization' weight). 

Another possibility is to replace ^/C^a by ^/C^a- Although this does not immediately do anything about the 

quadratic terms, if we also convert a portion vr of into obtain 



GrA =arA + (l-/u)^±±/rA + /-« 



i«rA, 



dz 



- + (M-^) 



+(1 - 2e)/C/CrA + 2^/2 



2TTri{fi-ri) 
1-4?? 



vr 



/CVrA + 



l-4r? 



rrA 



/C o /C/rA, 



(73) 



which requires C = 5 1 whereupon the two remaining equations become identical: /i ■ 



has many solutions. A particularly neat one is to take 77 
momenta' and vr 
BEFE's read 

GrA = -^rA + By A 



1-477 



, which clearly 



1 so that the primary objects are 'gravitational 



1 so that only two normal derivative terms appear in the 'bulk' term. Then /i = ^ so the 



-^rA = Tta + 



3rxx - r 
10 



/rA , -BrA 



1 



7nA 



1 SpAs^AEy.^^ 



10 dz 



(74) 
(75) 



Of course, it would make sense to particularly investigate the difficult geometrical problem '-BrA = on some 
embedded hypersurface' for such i?rA = corresponding to no quadratic terms, since by the same arguments as 
above, such a model would be sure to give braneworld physics distinguishable from that hitherto studied. 



4.5 Two Further Comments about Building SMS-type Braneworlds 

First, so far we have talked in terms of the 3^rA = ^rA — A/rA split of the matter contribution to relate our work 
as clearly as possible to its predecessors in the literature. However, from the outset [H] it was pointed out that 
this split is not unique. On these grounds we would prefer to work with the unique trace-tracefree split in which 
all the A/rA contributes to the trace part. The (4, 0) version of this split is used in Sec 6. 

Second, given a fixed type of bulk energy-momentum such as the Tab = of NKK or the Tab = ^dAB, 
then establishing an embedding requires the existence of a suitable compensatory characterization of the bulk 
geometry. The GR line of thought would be to only take results within such schemes seriously if they are robust 
to the addition of bulk matter fields. Of course privileged choices of bulk could arise from further theoretical 
input. We argue below that the theoretical arguments behind some privileged choices in the literature for the 
bulk energy-momentum are not convincing enough to anchor strongly credible physical predictions. One would 
rather require rigorous and general theoretical input following directly from some fundamental theory such as 
string theory. 

The vacuum choice of bulk of the induced matter NKK approach is always possible 
given the premises of the CM result. Then TZab = is claimed to be a complete geometrization of matter |4U| . 
However, the CM result holds equally well for any other analytical functional form of the energy-momentum. 
Furthermore, this approach considers only 1-component ('induced') matter; counting degrees of freedom shows 
that it cannot be extended to many important cases of fundamental matter. Whereas allowing for more extra 
dimensions could improve similar situations j4Ul I74j . unification requires geometrization of the fundamental 
matter laws themselves, whilst this 'induced matter' approach only geometrizes solutions of the EFE's coupled 
to matter of unspecified field dependence. 

^''For clarity, compare KK theory proper, in which the electromagnetic potential and Maxwell's equations are geometrized, as the 
frz portion of the 5-metric and the Grz equations of the KK split respectively, along with a scalar field. NKK theory could geometrize 
a more general vector field than the electromagnetic potential in addition to the generalization of the KK scalar field that is all that 
is usually considered. Now, this vector field is no known vector field of nature for its field equations does not in general correspond 
to that of any known vector field (although clearly the field equation contains the inhomogeneous vacuum Maxwell equations as a 



The Tab = ^9ab choice of bulk (an example of which is pure AdS) is clearly also always possible given the 
premises of a particular case of the generalized CM result. However, the motivation we wish to discuss is the 
argument for pure AdS bulks from string theory. This is not generally justifiable since firstly, bulk gravitons 
are permitted so the bulk geometry would generically contain gravity waves. Secondly, bulk scalars ought to be 
permitted since they occur along with the graviton in the closed string spectrum The content of the closed 
string spectrum thus places interesting restrictions on bulk matter rather than completely abolishing it. From 
the perspective of 5-d GR, evidence for the stability of vacuum or AdS bulks (and of any resulting physical 
predictions) to the introduction of suitable bulk fields would constitute important necessary support for such 
models and their predictions. 

Finally note that use of arbitrary smooth bulk Tab does not affect the form of the junction conditions since 
only the thin matter sheet contribution to Tab enters these. The conclusion of our second point is that there is 
no good reason not to explore at least certain kinds of bulk matter in order to have a more general feel for how 
these thin matter sheet models behave |75] . 

4.6 Sideways York and Thin-Sandwich Schemes with Thin Matter Sheets 

First we explain how our criticism of (3, 1; 1) methods in general in Sec 2.4 are applicable to the case with 
thin matter sheets. As regards the "evolution" w.r.t z, the issue of causality holds regardless of whether thin 
matter sheets are present, and the issue of well-posedness not being known will become particularly relevant in 
the study of sufficiently general situations in which rough function spaces would become necessary to describe 
the evolution of thin matter sheets (see Sec 6.1). The Campbell-Magaard scheme is of limited use in models with 
thin matter sheets not only because analytic functions are undesirable (and definitely inapplicable to sufficiently 
general situations) but also because the junction condition imposes restrictions on /CpA which prevents these 
being subdivided into the knowns and unknowns of Magaard's method. These restrictions remain even if one 
considers suitably well-behaved thick matter sheet models. In the next paragraph we outline how much we expect 
can be achieved with the (3, 1) version of the York method applied on the thin matter sheet, but recall that 
even if this does provide data sets, one is next confronted with the difficulties of the (3, 1; 1) "evolution" scheme. 
The main message is that the choice of methods which properly respect the difference between space and time is 
absolutely crucial. Thus, although at the simplest level (3, 1; 1) methods which build higher-dimensional bulks 
about the privileged (3, 1) worlds may look tempting, general attempts at proceeding thus are hampered by the 
required mathematical tools simply not being available and by these attempts not being in accord with the usual 
notion of causality. 

To date (4, 1) worlds built from (3, 1) ones have relied on very simple specialized ansatze, such as 

A) z-symmetric surfaces /CpA = with known metric /rAi whereupon the vacuum Codazzi equation is automat- 
ically satisfied and then R = is required from the Gauss constraint. 

B) /CrA = CfrA with known /rA; for example to obtain the Randall-Sundrum bulk jjj or slightly more general 
solutions 76 . Now the maximal subcase /C = of the CMC condition is a generalization of A), whereas the 
full CMC condition itself is a generalization of B). Moreover, now the metric is to be treated as only known 
up to scale. Thus one would generally only know the full metric of each model's (3, 1) world once the 'wave 
Lichnerowicz equation' for the embedding of this world into the (4, 1) world is solved. So one loses the hold 
from the outset on whether each model will turn out to contain an interesting (3, 1) world. Nevertheless, some 
of the (3, 1) worlds will turn out to be of interest. Furthermore, one should question the sensibleness of any 
ideas involving the prescription of full (3, 1) metrics if the most general technique available fails to respect such a 
prescription. This point is more significant for (3, 1) data than for (3, 0) data because conformally-related metrics 
have different non-null geodesies. For (3, 0) data no physical significance is attached to spatial geodesies, but for 
(3, 1) there are timelike geodesies which are physically interpreted as paths of free motion of massive particles. 
So a (M, hab) spacetime which is conformally related to (M, hab) is different physically (for example one could 
violate energy conditions the other one does not violate). One can get out of this problem by either attaching no 
physical significance to one's inspired guesses for (M, hab) or by hoping for unobservably tiny deviations between 
the geodesic curves of the two geometries. 

In the specific case of thin matter sheets, by the j.c's, A) implies that 3^rA = 0, whilst B) implies that y is 
a hypersurface constant. The maximal condition implies that y = whilst the CMC one implies that 3^ is a 
hypersurface constant. So whereas the maximal and CMC ansatze are more general than A) and B) respectively, 

subcase since KK theory would be included in this way within NKK theory). For sure, this vector field (or even more obviously 
the scalar field) is not capable of being a simultaneous geometrization of the individual vector fields of the standard model. If one 
were to treat these fields collectively, one throws away all that is gained by keeping their classical identities separate: the results of 
Weinberg-Salam theory and QCD that are obtained then by quantization. 



they are nevertheless restricted in this braneworld apphcation. Notice, however, that a number of interesting 
cases are included: vacuum, radiative matter and electromagnetic matter are all among the y = spacetimes. 

It is important to note that unlike the usual GR application, the choice of a hypersurface to be a brane is 
not a choice of slicing because localized energy-momentum is to be pinned on it. Almost all reslicings would fail 
to isolate this energy-momentum on a single slice. We know of no good reason why the brane should be CMC 
nor what value the CMC should take on it. However, at least this is a well-defined notion and it is substantially 
simpler to solve for than in general because of the decoupling of the Gauss and Codazzi constraints. 

As regards possible use of either of the two forms of the thin sandwich conjecture to treat branes, we first 
distinguish between thin sandwiches between 2 nearby branes and thin sandwiches which have a brane on one 
side and an undistinguished hypersurface on the other. One should be aware that non-intersection requirements 
may be different in these two cases, and also different from that of the original TSC setting of 2 unprivileged 
spacelike hypersurfaces. Second, one would have to take J^rA as unknown until it can be deduced from the /CpA 
evaluated from the TS procedure. Finally we caution that TSC schemes need not always exist. They require the 
"lapse" to be algebraically-eliminable, which for example is not the case for the analogous $ of the KK split. 

5 Applying (r, 1; 1) Methods to Remove Singularities 

We have argued in Sees 2-3 in favour of (4, 0; -1) schemes over (3, 1; 1) schemes. However, it is often suggested 
that (3, 1; 1) schemes may be useful since in certain particular examples (3, 1) spacetime singularities are 
'removed' by embedding into nonsingular (4, 1) spacetimes [SSI IHEj' Whereas embeddings are undoubtedly 
useful tools in the study of (3, 1) spacetimes (for the purpose of algebraic classification [771 |7S]), we ask what is 
the status of their use to remove cosmological singularities. Is it mathematically rigorous, generally applicable 
and physically meaningful? 

By 'singular' we mean geodesically-incomplete jTH]. In the study of singularities, timelike and null geodesies 
(t.g's and n.g's) are the only objects of fundamental importance since they are the curves privileged by the free 
travel of massive and massless particles respectively. The other objects associated with the study of singularities 
arise as a matter of convenience rather than being fundamental. 

Such objects include the expansion and shear of geodesic congruences, and curvature scalars. Because these 
objects are closely-linked to material in this paper, we focus first on them (Sec 5.1), explaining how they arise 
in the study of singularities and then demonstrating how the embedding and embedded versions of these objects 
can be very different (Sec 5.2). We then discuss the more important issue of geodesic incompleteness in the 
context of embeddings (Sec 5.3). An example is used to illustrate some of these points and others (Sec 5.4). A 
lack of rigour in singularity-removal claims is uncovered in Sec 5.5. We finally treat the case with thin matter 
sheets in Sec 5.6. 

5.1 Secondary Objects in the study of Singularities 

For (g, l)-dimensional cosmology, smooth congruences of past-directed normal t.g's with normalized tangents 
denoted by are considered. The decomposition 

Q 

Bafe = DaUb = -hab + (Tab + ^ah (76) 

q 

provides the expansion 6 = h'^^'Bab and shear aab = ^Jab)- twist uJab = ^[ab\ is zero for the normal congruences 
considered here, in which case Bab is an extrinsic curvature Oab- The corresponding normal Raychaudhuri equation 
[c.f ((221)] can now be considered. Although for > this would usually mean that a caustic develops, under 
certain global conditions a contradiction about the existence of conjugate points arises. Singularity theorems are 
thus obtained; in a cosmological context the simplest^^ of these is (Hawking, theorem 1 of |80j ) . 
For globally-hyperbolic (3, 1) GR spacetimes obeying the SEC and such that = C > on some smooth 
(spacelike) Cauchy surface S, then no past-directed timelike curve from S can have greater length than 

■^^ There are a number of other cosmological singularity theorems |14lllti| not only because some have weaker assumptions, but also 
because one wants to be able to treat a number of pathological cases. One such is the Milne universe: although this naively looks like a 
cosmology with its focusing of geodesies normal to t =const as t — > 0, it is merely a region of Minkowski spacetime. So one would not 
want to include the Milne universe among the singular spacetimes. The way out of this is to demand that physically-meaningful cos- 
mologies contain matter. This makes sense because the physical problem that may be associated with geodesies focusing is the possible 
pile-up of matter travelling along these geodesies leading to infinite densities. One succeeds in not including the Milne universe by use 
of singularity theorems hinging on the 'timelike genericity condition' Rabcdu'^u'' 7^ at least one point along each timelike geodesic, 
since the Milne universe is Riemann-flat. 



By the definition of Cauchy surface, all past-directed t.g's are among these curves and are thus incomplete, so 
the spacetime is singular. Such theorems generalize to {q, 1) spacetimes (provided that these obey analogous 
energy conditions) in the obvious way since Raychaudhuri's equation clearly behaves in the same way for all q 
and the global part of the arguments uses topological space methods that do not care about dimension. 

It must be noted that the singularity theorems are about existence whilst saying nothing about the nature of 
the singularity. Ellis and Schmidt began to classify singularities according to their properties [73], a difficult study 
which may never be completed [IH]- Below we consider only the most elementary type of genuine singularities: 
curvature singularities, for which at least one spacetime curvature scalar such as TZ or TZabR^^ blows up. 

5.2 Relating the Embedded and Embedding Secondary Objects 

Here we explain how knowledge of Gauss' hypersurface geometry renders it unsurprising that there are instances 
of singular spacetimes being 'embeddable' in nonsingular ones. This is because the behaviour of higher and lower 
dimensional objects used in the study of singularities clearly need not be related. As a first example, consider 
expansion. In the normal case, 9af3 = Qab = ~2a^3?^ from which follows 

e = ^e-eT^. (77) 

n 

So a blowup in 9 need not imply a blowup in 0. Thus the (4, 1) spacetime perspective on focusing of geodesies 
can be completely different from the perspective on some (3, 1) hypersurface. Figure 6 shows how one's (3, 1) 
notion of expansion generally corresponds to expansion and shear from the perspective of an embedding (4, 1) 
spacetime. In particular, what look like caustics or singularities in the (3, 1) spacetime could well correspond to 
no such feature in some surrounding (4, 1) spacetime. 

As a second example, consider curvature scalars, in particular the Ricci scalar. Preliminarily, Gauss' out- 
standing theorem is a clear indication that extrinsic curvature can compensate for differences between higher and 
lower-dimensional intrinsic curvature. Our simple idea is to consider the implications for rigorous embedding 
mathematics of the generalization of this to the case where lower-dimensional curvature scalars become infinite. 
From the form H19() for the (3, 1; 1) Gauss constraint, 

2p = 2n±± -n=^IC^ - IC^ olC'^ - R (78) 

clearly - at least for some (3, 1) worlds which have Ricci-scalar curvature singularities R — > oo - there will be 
surrounding (4, 1) worlds in which TZ (and the 5-d p) remain finite. For it could be that iC^ olC^ — |/C^ — > — oo 
suffices to compensate for R — > oo. Thus it could be shear, expansion or both that dominate the compensation 
for R — > oo. If this involves |/C| — > oo, it means a (4, 0) caustic forms in the (4, 1) spacetime corresponding 
to the (3, 1) spacetime's singularity. Yet this is not the only nonsingular possibility; if R — > — oo it may be the 
shear that compensates, including cases in which the blowup is pure shear. 

Finally, clearly the higher-dimensional singularity theorems hold so even if one were to succeed in excising 
singularities from a lower-dimensional model, one would typically expect singularities to occur elsewhere in the 
resultant higher-dimensional models. The nature of these higher-dimensional singularities may not be the same 
as that of the excised lower-dimensional ones. So, interestingly, by allowing extra dimensions, one would be even 
less certain of the character of singularities. 

5.3 Embeddings and Geodesic Incompleteness 

The fundamental importance of geodesies is problematic for embeddings since in general the (3, l)-dimensional 
geodesies are not included among those of an embedding (4, l)-dimensional spacetime. There is then the dilemma 
of which of these sets of curves should be the physically privileged ones and thus be the set of curves whose 
extendibility is in question. 

If one wishes to postulate (4, l)-dimensional GR, in addition to the (4, l)-dimensional EFE's one must require 
that the matter follows the (4, l)-dimensional geodesies. Then the issue of extendibility of the original (3, 1) 
geodesies becomes irrelevant since they are unprivileged curves in the (4, 1) spacetime. But this is at the price of 
altering what the effective (3, l)-dimensional physics is: often it will be of affine-metric type and thus in principle 
distinguishable from ordinary (3, 1) physics. One cannot then necessarily claim that the recovered (3, 1) physics 
is as close to the usual (3, 1) physics as is widely claimed in the braneworld literature 

If these do not include the (3, 1) geodesies, then the (3, 1) question of what happens to these becomes 
physically irrelevant. One would not then be extending one's (3, 1) geodesies but replacing them with distinct 



(4, l)-dimensional geodesies, whieh might be viewed as unsatisfactory. There is the further comphcation in 
braneworld models that one might wish for (3, 1) and (4, 1) geodesies to both be relevant, as the paths of 
brane-bound and brane-unbound matter or gravity waves. It is then not clear what is meant by a singularity - 
exactly which curves are supposed to be incomplete? 

We end with two loose possibilities about incompleteness which might sometimes be related: the (3, 1) 
singularities might be interpretable in terms of a (3, 1) hypersurface that becomes tangent to the characteristics 
of the bulk equations, and incompleteness might correspond to t.g's and n.g's being forced to exit the (3, 1) 
hypersurface. The second possibility is a variant on the theme of (4, 1) and (3, 1) geodesies whereby the (3, 1) 
geodesies are pieces of (4, 1) geodesies which are extendible only by replacing a piece of the original foliation 
with a new one that extends into what was originally regarded as the extra dimension. 



5.4 Embedding Flat FLRW in Minkowski 

The following example illustrates many of the above points and leads us to further comments. Consider a (4, 1) 
spacetime with metric (see e.g |32j ) 



QAB = diag(5tt, g^^) = diag \-z , t^z^-id^p, _ 1 , q>l. 



(79) 



This is simple to treat because: 

1) If we foliate it with constant z hyper surf aces, a portion of each z = const hypersurface has induced on it a flat 
FLRW cosmology. In particular, the coordinates have been chosen so that the z = 1 hypersurface is the FLRW 
cosmology with scale factor ti . We restrict attention to g < 3 due to the DEC. For q 2 there is a (3, 1) Ricci 
scalar curvature singularity (g = 2 is the radiation universe, whence R = S = 0). 

2) The (4, 1) spacetime is in fact Minkowski. 

First, foliating the flat FLRW cosmology with constant t surfaces, ^ = 4- — > cxd as t — > 0, so there is 

gt 

(3, 0; -1) focusing as the Big Bang is approached. Furthermore, only focusing occurs: = 0. Next foliate 
the (4, 1) spacetime with constant t surfaces. The (4, 1) spacetime contains many FLRW worlds on constant-z 
surfaces. Build a congruence by collecting the FLRW geodesies on each z = const slice. Then the (4, 0; ~1) 
expansion Q = also blows up as t — > oo, but there is also a blowup of the corresponding (4, 0; -1) shear: 

so we have case 2 of Fig 7. Also, whilst still both — > oo and B — > oo in this case, we know that the former 
viewed from within the z = 1 (3, 1) hypersurface corresponds to a genuine (3, 1) singularity, whilst the latter in 
(4, 1) would be a mere caustic. 

Second, consider the z = 1 slice of the (4, 1) spacetime. Here the (3, 1; 1) expansion and shear are 

/C = ^ , /C?A = -^diag(3, tU^s). (81) 

Thus for the physical range, both blow up as t — > 0. So the spacetime includes a z = 0, t = point, at which 
there is a (4, 0) caustic. These blow-ups combine in /C"^ o /C"^ — |/C^ to cancel R for all values of q; for q = 2 the 
shear and expansion contributions exactly cancel each other. The unphysical q = 4: case is an example of a pure 
shear blowup. 

Although the congruence two paragraphs up is built to naturally include the geodesies of all the included 
FLRW universes, these turn out not to be the (4, 1) geodesies (nor pieces of them). Rather, by considering the 
coordinate transformation between (|79j) and standard (4, l)-dimensional Minkowsi coordinates, it is easy to show 
that the Minkowski geodesies pierce the z = constant surfaces that are the FLRW universes (more significantly 
in the early universe than today). As t — > the z = constant surfaces approach the null cone so the example 
is an illustration of the foliating (3, 1) hypersurfaces becoming tangent to the characteristics at the point of 
interest. Also the foliation breaks down as t — > because the family of hypersurfaces of constant z intersect at 
t = 0. All these points are illustrated in figure 6 b). 

Finally, the example cannot be taken to be typical by the (4, 1) version of genericity. Indeed, we know that not 
many (3, 1) spacetimes can be even locally embedded into (4, 1) Minkowski fTHl. The value of particular examples 
is limited to exhibiting possible behaviours. One would require this to be upgraded to the study of large classes 
of examples to assess probable behaviours. Indeed, one would expect higher-dimensional cosmological models 



to exhibit a range of behaviours distinct from that of the Big Bang (including nonsingularness). This fohows 
from knowledge of particular examples (from the algebraic study) of inhomogeneous cosmologies in standard 
GR We especially note that if highly-symmetric cosmologies are haphazardly embedded, one would often 
expect the resulting higher-dimensional models to be less symmetric i.e more anisotropic and inhomogeneous 
than the original lower-dimensional models [as naively suggested by the appearance of (4, 0; -1) shear from 
models containing no (3, 0; -1) shear]. So the resulting higher-dimensional models are expected occasionally 
to exhibit more unusual behaviour than the lower-dimensional models from which they arose, and this includes 
the possibility of nonsingular cosmologies. Indeed, we would expect that an increase of dimension increases the 
variety of possible behaviours in inhomogeneous cosmology. But the bottom line in the study of 4-d and 5-d 
singularities is the same: one is ultimately interested in generic behaviour and this requires more sophisticated 
means of study than building individual particular examples. 

5.5 Nonrigorousness of Singularity Removal by Embedding 

There are arguments against the use of embedding results toward making general statements about singularity 
removal. Our Sees 2.4 and 3 provide arguments against the use of the CM result (some of which also hold 
against the use of further embedding theorems). In particular, given a spacetime, there are so many possible 
embeddings. Were nature higher-dimensional, why should it choose a particular nonsingular embedding out of 
an infinity of (nonsingular or singular) embeddings that mathematically exist? The more extra dimensions are 
present, the more severe this nonuniqueness is. This makes singularity removal by such embeddings physically- 
questionable. The CM scheme, and other embedding schemes such as the 10-d Minkowski embedding theorem 
j77j ) . are local. This localness includes the notion that the theorem is only applicable to a small region of the 
original manifold. We question whether this need include the approach to a singularity, since these are edges of 
the lower-dimensional spacetime and thus have no neighbourhood. Singularities are global features of spacetime. 
However even some of the global embedding theorems (which require very many extra dimensions) are stated 
for the analytic functions 77 . But low differentiability may be typical in the approach to singularities [IS]- It 
becomes much harder to make any general statements once one accepts that spacetime is not analytic! 

One may view at least some (3, 1) singularities as projective effects due to taking a badly-behaved foliation 
(for example a foliation that becomes tangent to the characteristics), manifested by the formation of caustics or 
pure shear blow-ups. The question is why there is confusion in reversing this projection procedure to embed 
(3, 1) singularities into nonsingular (4, 1) spacetimes. If an embedding of the (r, s; e) type such as those 
considered in this paper is to be used, it is mandatory that the region of hypersurface to be embedded be 
entirely of one signature s (the CK theorem demands a nowhere- characteristic lower-dimensional manifold). 
Also, from the lower-dimensional perspective, the singularity cannot be included in the set of points on which 
data is prescribed because it is an edge of spacetime, not a point of spacetime. Were one to try and include it 
by providing data 'right up' to that edge, one would expect that the data would become badly-behaved e.g some 
components of /CpA could be forced to be infinite. Under the various circumstances above, embedding theorems 
become inapplicable to reconstruct the higher-dimensional spacetime. 

5.6 (n, 1; 1) Singularity Removal and Thin Matter Sheets 

Finally, we have two comments about the application to thin matter sheets. First, there is less scope for singularity 
removal along the lines of (|78|) . For, as it is the energy- momentum on the brane, 3^rA is presumably finite for 
a nonsingular (4, 1) world so /CrA is finite and so cannot cause blowups in fC^ o KL^ — |/C^. However, blowups 
in this quantity can still occur if the inverse metric is badly-behaved (corresponding to / = 0). Note however 
that the application of the CK theorem to GR requires this not to be the case everywhere within the region of 
applicability. 

Second, one would of course require the investigation of more elaborate (e.g 'warped') embeddings than the 
example discussed above in order to investigate whether some embeddings lead to 5-d geodesies that exit the 
brane. If the (3, l)-dimensional geodesies on the brane are not included among those of the (4, 1) dimensional 
bulk, it is not clear at all what is meant by 'singular' since one then one would have to simultaneously consider 
some notion of extendibility two congruences of privileged curves of different dimensionality. 

Under the decomposition the 5-d geodesic equation + T^bcx^x^ = becomes 



+ r^Asi^i^ + 2/C^A = 



(82) 



(83) 



if one upholds the use of normal coordinates. In order for ()82() to reduce in all cases to the (3, 1) geodesic equation 
on the brane, one requires everywhere on the brane X"*" = or /C^a = 0. By ()83() . maintenance of i;-*" = along 
all geodesies is impossible unless /Cas = 0. Thus in all cases /Cae = is required on the brane. But this means 
that /Cas^ = ^AE~ or equivalently 3^as = i.e the absence of a brane. Thus in SMS-type braneworlds the 
(3, 1) geodesies are not generally among the (4, 1) geodesies. So if one wishes to adhere to the postulate that 
bulk-traversing fields follow (4, 1) geodesies and fields confined to the brane follow (3, 1) geodesies, one should 
accept that the restriction of the former to the brane generally does not coincide with the latter. Thus there may 
be observable consequences in the (3, 1) physics: freely-falling bodies composed of confined matter would follow 
distinct paths from freely-falling bodies composed of matter capable of traversing the bulk but which happens 
to be moving along within the brane, which represents a violation of the equivalence principle, and the causal 
structures within the brane due to confined photons and due to gravitons which happen to be moving along 
within the brane would generally be distinct. Also the above new conceptual difficulty about the definition of 
singularities is thus indeed relevant to brane cosmology. 

6 Application of (n, 0; —1) Methods to Thin Matter Sheets 
6.1 Difficulties and Modelling Assumptions 

We are interested in the study of (n + l)-dimensional GR models with thin matter sheets using only well- 
studied, well-behaved mathematical techniques which also make good physical sense and are applicable to full 
EFE systems. In our view a good way of achieving this is via the (n, 0) IVP, followed where possible by the 
heavily-protected (n, 0; -1) CP. 

We provide a hierarchy of modelling assumptions together with associated sources of difficulty, and argue 
that evolution is more problematic than data construction. We concentrate on the latter, making use of all the 
attractive features of the York method mentioned in Sec 3, albeit this elliptic mathematics is in a new setting 
(see Sec. 6.2-3). The differences come about due to the following difficulties of relevance to the hierarchy of 
modelling assumptions. Firstly, the roughness of function spaces required in the rigorous treatment of certain 
thin matter sheet models. Secondly, the requirement of novel asymptotics and boundary conditions (b.c's). 

Note that the increase in dimensionality itself does not affect much the data construction nor the study of 
well-posedness and stability of the GR CP, although there is some dimensional-dependence when the rougher 
function spaces come into use, due to the following two contrasting results. 

The first result is that in all dimensions, the presence of thin matter sheets means that the metric is rougher 
than (the functions with continuous second derivatives) because of the jump in the Riemann tensor at the 
sheet. This corresponds to the metric being too rough to belong to the Sobolev class H"^. 

The second result is that the definition of a product^^ for Sobolev spaces for the EFE's in harmonic coordinates 
according to Hughes, Kato and Marsden (HKM) |47[ I49j requires the Sobolev class of the metric to be no rougher 
than H^^^ [i.e H'^ for standard GR and for (4, l)-dimensional GR]. Thus HKM's mathematics, the strongest 
used for both the GR IVP and CP in the review ^J, is not generally powerful enough to deal with thin matter 
sheets. 

More powerful mathematics is thus desirable, and more may well appear over the next few years (see the 
program starting with HE])- Realizations of part of the hierarchy of models below would be an interesting 
application for this sort of mathematics, which further motivates the study of increasingly rough Sobolev spaces, 
since the higher the dimension the greater the improvement to the established HKM mathematics is required 
for a full study of thin matter sheets. We must also mention that in spaces rougher than there is no 
guarantee of geodesic uniqueness ^| , so complete conventional physical sensibleness is generally lost no matter 
what rougher Sobolev class results are proved. This shortcoming is clearly dimension-independent (from the 
dimension-independence of the form of the geodesic equation). At least one need not go below whereupon 
geodesies need not even exist! 

The asymptotics problem 82 follows from supposing that 'AdS' bulks are desirable for string-inspired sce- 
narios. The study of 'AdS' bulks by relativists would entail the study of asymptotically AdS bulks, to permit 
a more general study of disturbances close to the thin matter sheet. The following comments are in order here. 
First, in contrast to dS, AdS is infinite and hence possesses asymptotics. Second, whereas we will for the moment 
seek to avoid the function space difficulty, the asymptotics problem will remain in the 'tractable scenarios'. In 
the usual GR, the use of AdS asymptotics for the application of the York method to small-scale astrophysics was 
neglected because of the negligible effect of a cosmological constant on such small scales and because of DOD 

■^'^This is a necessary complication because the EFE's are nonlinear. 



arguments. The braneworld application is however substantiahy different, so one might in this case have to study 
the York method with 'AdS asymptotics'. This could affect the tractability of the Lichnerowicz equation, and 
also whether the existence of the crucial maximal or CMC slices is as commonplace for these new spacetimes 
as it was for the (3, 1) GR ones. Furthermore, since one is interested in infinite planar branes, the asymptotics 
is directional, 'perpendicular' to the brane. However, in general there is no such notion of perpendicularity in 
GR. Whereas one can locally define geodesies and draw hypersurfaces perpendicular to them, this procedure in 
general eventually breaks down, for example due to (spatial) caustic formation. 

We start our hierarchy by envisaging the most general situation possible for thin matter sheets within the 
framework of GR. It is entirely legitimate to construct a (4, 0) initial data snapshot with whatever shape of thin 
matter sheet, but it is not legitimate to assume that any features of this are maintained over time unless it can 
be shown that the full evolution equations robustly maintain these features. Strictly, this is an initial-boundary 
problem, such as occurs for water waves or for the surface of a star. This is a very hard and quite new problem 
in GR j45| i84j . Thus the evolution step is particularly hard both in its full generality and in the justifiability of 
simplifications such as a fixed boundary. 

Consider a thin matter sheet in an asymmetric (not Z2 symmetric) bulk. One would ideally want to follow 
objects that crash into and possibly disrupt the thin matter sheet. One can then imagine asymmetric crashes 
which might punch though the thin matter sheet. The sheet could thicken or disperse with time. It could emit a 
significant amount of gravitational radiation, which could moreover include gravitational shock waves that spread 
out the character that initially pertained only to the thin matter sheet and not to the smooth surrounding 
bulk. That is, the character of a thin matter sheet data set J could typically infect the whole of its causal 
future J~^{J) (Fig. 7). In simple words, why should the junction remain unscathed? Why should the bulk in 
the immediate vicinity of the junction at later times be simple, smooth or known? Addressing these questions 
goes beyond the reach of present-day techniques. Rather we note that some usually-tacit assumptions are used 
to definitely avoid these difficulties, and that we do not know whether any of these are justifiable in the physics 
governed by the GR evolution equations. One can presuppose a thin matter sheet exists at all times to prevent 
it being created or destroyed. One can allow only bulks with regular thin matter sheet boundaries to preclude 
shock waves. One could also preclude asymmetric crashes by presupposing Z2 symmetry. 

Even if discontinuous emissions are precluded, one can imagine smooth emissions and absorptions (symmetric 
or asymmetric) whereby the bulk interacts with the thin matter sheet. This would entail material leaking off 
or onto each side of the thin matter sheet. This can be precluded by the presupposition that the thin matter 
sheet's energy- momentum resides at all times on the sheet, which is encapsulated by the well-known 'equation of 
motion' of the thin matter sheet |12l I19j . This is a significant restriction on the dynamics of thin matter sheets 
in GR, often carried out in the name of tractability but which may not be realistic. 

In the Randall-Sundrum scenario, the brane can be placed anywhere in the bulk without affecting it. We 
strongly suspect that this is not a typical feature in GR-based braneworlds, nor is it desirable since it dangerously 
marginalizes the ontological status of the bulk. We take this feature to be too specialized to be included in the 
developments below. 

We close this section by briefly discussing modelling assumptions which avoid some of the above difficulties. 
We then implement these in two particular classes of tractable problems. The desirability of each modelling 
assumption strategy and of other modelling assumptions from a string-theoretic perspective is discussed in Sec 
7. 

I) The choice of Z"^ symmetry may well look and indeed be arbitrary from a GR point of view and indeed 
from a string-theoretic point of view (see Sec 7). But without such a choice one simply cannot establish any 
specific embeddings if there are thin matter sheets (in the GR study of stars one requires the absence of surface 
layers to perform matchings). Thus one studies the restricted case with Z"^ symmetry, in which one integrates 
only up to the b.c provided by the j.c. In the IVP part one can assume whatever configuration of thin matter 
sheets. Thus this step is more justifiable than the subsequent evolution, which may require some of the above 
ad hoc assumptions as to the existence and good behaviour of the thin matter sheet at all future times. So for 
the moment we are merely after the construction of (4, 0) initial data up to the junction. By themselves these 
data are useful in addressing issues such as the shape of the extensions of black hole horizons into the bulk, but 
the staticity and stability of these configurations remain unaddressed in the as-yet intractable evolutionary step 
of the problem ^] . 

II) One could choose instead to work with thick matter sheets i.e ones that are finitely rather than infinitely 
thin [Mj. Then the above troubles with the function spaces need not occur, since the walls would then not be rough 
in the above sense, although they could still have the tractable level of roughness that is able to accommodate 
astrophysical objects just as in ordinary GR. It is then easier to envisage the study of the evolutionary step of 



the thick matter sheet problem quite high up our hierarchy, since this more closely resembles the stellar surface 
problem. 

Ill) One could investigate closed shells rather than open sheets in order not to require directional asymptotics. 



6.2 More Techniques for the (n, 0) IVP 

We cast the initial value constraints as 

A<j) - Nonlin(0) = , Nonlin = -c(l){R - Mcp" + /i^'' - Ipcj)") , M = o (84) 

DaQ'^^' = f (85) 
4 1 J arbitrary constant = 



n-2 ' {l-n)b ' 1 ^ //O. 



(86) 



We first consider the j'^ = case. Then by definition Q^"-^ = 0TTafe g^g^^gg \yl order to solve the Codazzi 
constraint. Furthermore, in that case a particular scaling of is not required to keep the Codazzi constraint 
conformally-invariant, so preserving the DEC does not enforce a particular scaling of p. So we are left to solve 
just the Lichnerowicz equation with the scaling of p to be determined by our convenience. 

We begin by considering the linearized Lichnerowicz equation. Set (p = (po + ecpi, p = po + spi in H84|) and 
equate the 0{e) terms to obtain 



a+a 



[A - Lin(0o)]0i = lnh(0o) , Lin(0o) = -c[R - (1 + c)M4,l + p^{l + b)4 - 2(1 + a)p4>l] , Inh(,/.o) = 2cpi<; 

(87) 

If asymptotically hab — > h^f, and (po — > (po, one obtains [A^ — Lm{(f)Q)](f)i = Inh(^o") from which the asymptotic 
behaviour of cp for the full equation may be obtained. Eq. (|87() is a linear elliptical equation^^ 

[h'^didj - - Lm{cPo)]cPi = Inh((/<o) (88) 

(where F* = h^'^T'^jk). Many b.v.p's for such equations are well-studied O |HS1 IHEl- To give an idea of the 
techniques available, consider the Dirichlet problem: (|88p in S together with b.c (pi = f on dT,. On sufficiently- 
benign function spaces, the well-known uniqueness proof based on Green's identity holds, provided that Lin > 
everywhere. This uniqueness can also be proved from the maximum principle, the use of which applies to more 
general function spaces. One then proves existence from uniqueness |86J , by setting up a sequence of problems 
interpolating between the flat-space Laplace equation and (|88j) . One can then interpolate between quasilinear 
elliptic equations such as the full Lichnerowicz equation and their linearizations. Then one has topological 
theorems about the fixed points of suitable maps, such as Leray-Schauder degree theory [HEl, to prove existence 
of the solutions to the quasilinear elliptic equation's Dirichlet problem. Similar techniques are available for 
further b.v.p's |86] . 

Thus the chances of having a well behaved b.v.p for the Lichnerowicz equation are much enhanced if the 
corresponding Lin > 0. To this end, notice that M and p"^ are always positive whereas R is of variable sign. 
Contrary to what York assumes in the ordinary GR context, we must permit p < since our application requires 
a negative bulk cosmological constant, which further complicates the analysis of which cases are guaranteed to 
be well-behaved. There are three further simplifications one can try: 

i) Conformal flatness, which replaces A by the flat-space Ap and wipes out the arbitrary-sign R term. 

ii) Maximality: ^ = removes the highest-order term of Nonlin((/)). 

iii) Setting M = 0. If this is combined with ii), one has time symmetry Qij = 0, which is regarded as overly simple 
in the standard applications of numerical relativity as it corresponds to the absence of gravitational momentum 
on the slice S. 

One can sometimes also use a trick: to scale p like M, R or p'^ and then fix it so as to cancel with this term. 
Because we are interested in p < 0, only the first two of these possibilities are available. 

In the j"" 7^ case, the cancellation trick above is not available if one wishes the procedure to respect the 
DEC. Simplification i) is convenient in this case. Whereas it is important to develop this case, it lies beyond the 
scope of this paper. 

■^^Note however that h^-' and F* have geometrical significance: our linear elliptic equation is on an underlying n-manifold E (generally 
with boundary). We thus caution that if S cannot be deformed to flat space, the interpolation methods between the 5R" Laplacian 
on a portion of and our operator cannot apply. To avoid this trouble, one could study a sufficiently small region of E (as E is a 
manifold and hence locally or restrict the allowed topology of E to be the same as that of a portion of K". 



6.3 The (4, 0) IVP with Thin Matter Sheets 



We work with the unphysical hne element ds^ = dz'^ + e^dxadx'^ for some known trial function w = w{r,z), 
so as to model a S'^-symmetric object such as a black hole or a star on the brane. For the = case, the 
IVP essentially reduces to the solution of the Lichnerowicz equation for some (j) = We begin with the 

Tab = case of (the spatial projection of) ()5U() . We first find the b.c's for the Lichnerowicz equation, then provide 
algorithms and then comment on the underlying analysis. 

Note that the b.c's are to be imposed on the physical metric and then written in terms of the unphysical 
quantities to have b.c's in terms of the objects we start off with and the unknown conformal factor cj). We begin 
with the inner conditions and then discuss asymptotic outer conditions. To establish this sort of inner radial b.c. 







it is commonplace to use an isometry [S7]. For example, one could use hab{xa) = hab{—Xa) =^ 

leading component-by-component to the homogeneous Neumann b.c (the so-called reflection b.c 



dr 



r=0 



dr 



0, 



r=0 



and the restriction ^ = on the valid form of the known function. 

One might worry in the case of black holes that r = is singular. For ordinary GR, one excises the singularity 
by use of an inversion-in-the-sphere isometry about some throat at radius a of the black hole within the apparent 
horizon (See Fig 8a). This gives an inner Robin condition |87j 



— + — (; 
dr 2a 



0. 



(90) 



For braneworld black holes one does not know how far the black hole extends into the bulk. The main point of 
the study is to find this out (pancakes versus cigars). Our idea is then to guess a profile r = /(z) along which to 
excise. About each point on r = f{z) an isometry in the corresponding 2-sphere may be applied (fig 8b), leading 
to an inner Robin condition 

1 



dr + 2f{z) 



= 0. (91) 

r=f{z) 

The profile could be chosen so that it matches up smoothly with the r = and z = 1 boundaries. Once the 
problem is solved, one can find out whether this profile was a good choice or not (see fig 8c). By picking a 
1-parameter family of r = f(z) curves, one could then iterate until a satisfactory profile is found. 
For the inner-z (z = 1) b.c, we impose the j.c's ()41() and H49|) . which now read 



K 



af3 



Y- 



a/3 



lap 



(92) 



(N.B these are imposed on the physical quantities). Imposing Z2 symmetry, the second of these becomes 




4\ 



(93) 



z=l 



by use of the definition of the physical K^js in normal coordinates and then factorizing out the Cap- This Z2 
isometry b.c is more complicated than usual because of the presence of the thin matter sheet. 

That the Z2 symmetry is termed a refiection appears to generate a certain amount of confusion. This is 
because there is also the refiected wave notion of reflection, taken to give a Neumann b.c (such as in the r = 
condition above). In the case of Shiromizu and Shibata (SS) [21] (see also 1251). who furthermore treat our known 
w as their unknown, and our unknown (p as the perfect AdS (/> oc ^, it has the nice feature that the other 2 terms 
in H93|l cancel leaving a Neumann b.c 

= 0, (94) 

=1 



z+q{r) 



which 



dz 

which might look like a reflection in both of the above senses. However, this is only possible for (f) oc 
may be somewhat restrictive. Note also that (jHU is an unusual notion of reflection in this context, since it is 

■^**One has no choice but to work with apparent horizons rather than event horizons if one is given a single spatial slice. At least 
the apparent horizon lies within the event horizon in GR. 



of the conformally-untransformed (3, 0)-metric, whereas the natural notion of reflection isometry would be for 
the conformally-transformed (4, 0)-metric. Thus we conclude that the presence of the notion of reflection H94|) in 
SS's work is happenstance: it is separate from the Zi notion of reflection encapsulated in (|92() and no analogue 
of it need exist in more general situations. 

Our b.c. (|93jl contains a reflection part the other part being a (nonlinear) absorption. This has the 
following implications. First, the presence of absorption terms is interesting since pure reflection is the boundary 
condition of a perfect insulator in the potential theory of heat. This suggests that pure reflection b.c's such 
as for pure AdS have built-in non-interaction of the bulk with its bounding brane, whereas our b.c may lead 
to (evolution models) with brane-bulk interactions. Whereas neither the planar symmetry in (4, 1) spacetime 
of the underlying thin matter sheet nor the 5^ symmetry within the (3, 1) thin matter sheet of a simple 
astrophysical object generate gravity waves, we expect that spherical bumps on approximately planar branes 
in (4, 1) spacetimes to admit gravity waves, potentially giving rise to instabilities in the approximately AdS 
bulk braneworlds in which the branes contain astrophysical objects. Our b.c could permit the modelling of such 
interactions. Second the nonlinear absorption (much as arises in the theory of heat with temperature-dependent 
conductivity) is likely to complicate both analytic and numerical treatments of our b.v.p. 

The first point above is one reason to favour our approach over that of SS (subject to overcoming the 
complications due to the second point). We now provide a further reason why the approach we present should 
be favoured on the long run. Whereas one could try to generalize SS's work by keeping the notion that w be 
unknown and (\) known but not whereupon our b.c (|93() is interpreted as an inhomogeneous Neumann b.c in 
tD, the main trouble with SS's work is that there is no good reason for their method to be directly generalizable 
away from its restrictive assumption Ojj = 0. Whereas treating w not <^ as the unknown may give nicer b.c's, it 
is not tied to a method known to be amenable to less trivial solutions of the Codazzi equation than O^- = 0. Our 
proposed method seeks a workable extension to the significantly more general case 0^- = ©ij"'", and speculatively 
to completely general ©jj. 

Another nice feature in SS's work is that the matter is scaled so that the linearized p.d.e obtained is a Poisson 
equation. This p.d.e is simply invertible to obtain detailed asymptotics |HH1I23- For our scheme, the linearized 
equation ()87() is of Helmholtz-type rather than Poisson, complicating such a procedure. We merely demand 
instead that — > ^ for z large wherever possible so that our models are 'directionally asymptotically AdS'. 
As for large r, at least on the brane, one can impose asymptotic flatness as r — > oo. It is less clear what one 
should prescribe in this respect off the brane. A more detailed study of the asymptotics for our b.v.p's should 
be required as part of the actual construction of examples of data sets.^^ 

Our proposed b.v.p is thus the mixed, nonlinear problem 

= Nonlin(0) , 

l^ + f (^-^)=Ofor. = landr>/(l), <^ ^ i as . oo , ^^^^ 

1^ = for r = and z > /~^(0) , (j) — ^ 1 as r — > oo 
. § + ^J = Oonr = fiz). 

One can use simplifications i) and/or ii) on this and still have a more general case than SS. Simplification i) 
might be plausible because 5-d AdS clearly admits conformally-flat spatial sections [ds'^ = -^{I'^dz'^ + dx^dx"')] 
or maybe not, since it is regarded with suspicion in the (3, 1) GR 2-body problem The trick at the end 
of Sec 6.2 can not be used to obtain particularly simple examples due to the following argument. Unless one 
puts p (X A (constant) and has this maintained by non-scaling (i.e a = 2), it is overwhelmingly probable that the 
emergent A is not constant. Thus our freedom in a is used up to ensure that A is constant. Our algorithm for 
solving this problem is as follows: 

i) Prescribe the following unphysical quantities: the metric hab and matter density p. 

ii) Pick any suitable O'^'^*-' to solve the Codazzi constraint. It follows that M is known. 

iii) Thus we can attempt to solve our b.v.p for the Lichnerowicz equation to obtain (p. 

iv) Then we can compute the physical bulk metric hab and induced brane metric e^/? of our snapshot. 

This assumes that we have correctly guessed the profile. One would now check whether this is the case by 
solving for the apparent horizon. If this is unsatisfactory (fig 8c) then one would repeat the algorithm with 
adjusted profile. 

For nonvacuum branes {Tab 7^ 0), a similar argument enforcing a = 2 holds. The j.c may now be split into a 
trace b.c part, 

■^^Particularly because numerical integration is done on large but finite grids, subleading order asymptotics are helpful. 



dz 2e\dz 6 
and a restriction on the tracefree part of the matter on the brane, which is 



= (96) 

2 = 1 



Yip = (97) 

for the metric ansatz and coordinate choice used. Our algorithm now becomes 

i) One now requires the prescription of hac everywhere and of Y . 

ii) Pick any suitable O'^'^*-' and hence M. 

iii) Attempt to solve our b.v.p for the Lichnerowicz equation. 

iv) Then we can compute hab and iap- In this particular case, the simultaneous imposition of normal coordinates 
and an isotropic line element forces the restriction (|97() on the braneworld matter. This is because one is 
applying too many coordinate restrictions; in this light the habitual practice of working in braneworlds using 
normal coordinates may be seen as a poor choice of gauge. If the above coordinate choices are not simultaneously 
made, the algorithm would contain further nontrivial equations [in place of (|97|) ] to solve before the braneworld 
matter content can be deduced. 

We assume that this nonvacuum application is for a star or clump of dust on the brane in which case there 
is no need to excise a corner with a profile as was done above to deal with the black hole singularity. 

Notice the lack of control of the physical metrics characteristic of theoretical numerical relativity. Our non- 
scaling of the matter at least gives visible control over the matter. The sensibleness of doing this is in fact tied 
to the negativity of the bulk p: for p > and a > — 1, Lin receives a negative contribution with its tendency to 
encourage ill-posedness, whereas for predominantly /c < 0, the danger lies in a < —1. 

A simpler example of b.v.p set up along the lines we suggest is that of Nakamura, Nakao and Mishima |25j . 
Their metric is cylindrical not spherical and they have a slightly simpler version of ()93() . Their case is still time- 
symmetric. The linearized equation is then precisely Helmholtz, facilitating its inversion and the consequent 
more detailed knowledge of the asymptotics. 

We also propose the 2-brane version of the above, in which the z — > oo condition is replaced by another 
'parallel' brane boundary at z = 1 + Z. This problem is close to that considered by Piran and Sorkin [55]. In 



contrast however, their study involves a periodical 
isometry-based reflection Neumann condition ^ 



ly-identified fifth dimension to which there does correspond an 
= 0. 

=1 

We can provide a local (in function space) uniqueness proof for our proposed type of b.v.p, at the same level 
as that for standard GR black hole data in |S7j. In the conformally-flat case, suppose that ipi = ip2 + u for u 
small. Then the homogeneous linearized b.v.p in u is applicable. We then have 

/ (I Vup + Uuu')d'x = fi\Vu\' + uAu)d'x = i updS = ^^f ''^%^JrV^%d<^ (9^) 

by Green's theorem. For u — > at least as fast as -, the first surface integral tends to zero. For u — > at least 
as fast as ^ the second surface integral tends to zero. The fourth integral is zero by the Neumann reflection b.c. 
Upon applying the linearization of the nonlinear b.c (|9!-{j) . 



du u f dw ~ 

Tz ^ 2 [j; - ^'^^'^ 



0, (99) 



the third surface integral is non-positive provided that the restriction ^ > holds. Thus since Lin > u must 
be zero. For the non-conformally flat case, this argument requires R> 0. The argument is not seriously changed 
in the black hole case when a corner is excised by a profile with a Robin condition on it. 

Stronger existence and uniqueness proofs for our proposed problems are complicated by four factors: 
unboundedness, (albeit of the most benign kind), the boundary having corners, mixed b.c's (i.e a piecewise 
prescription on the boundary) and a portion of these being nonlinear. The first is however of the simplest kind 
and the second and third are commonplace. The last is less usual but perhaps not so bad either because by 
A(p = iy^^^'' 'Bx'^ problem can be written in divergence form, for which b.v.p's with nonlinear b.c's are 

treated in Ch. 10.2 of 86^. The trouble is that the treatment there, unlike here, is neither for mixed b.c's nor for 
an unbounded region with corners. Thus, the strongest point we can make at present is that there is good hope 
of obtaining existence and uniqueness results even for quite rough function spaces by more-or-less conventional, 
entirely rigorous mathematics for our proposed problem. This should be contrasted with the hopeless state 
of affairs with sideways prescriptions! Once one begins to get good numerical results, it becomes worthwhile 



to explicitly prove the well-posedness of the method used, in order to support those results and the ongoing 
production of more such results. In this particular case, these results would serve to support and understand the 
crucial numerical step iii) in the above algorithms. 

To permit the even greater generality required to have momentum flows (j" 7^ 0), by our route one is forced 
to forfeit the DEC control over the matter. This is because at least the A part of p cannot scale in accord with 
the scaling of j"'. We leave this further development, which additionally requires posing and solving the b.v.p 
following from the Codazzi equation, for a future study. 



6.4 The (4, 0) IVP with Thick Matter Sheets 

Many of the difficulties with evolution discussed above stem from the thinness of the matter sheets causing 
function-space-related problems. But ordinary physics should not be sensitive to which function space is used. 
One would hope that models with thick matter sheets would be more amenable to study, be good approximations 
to the thin matter sheets and in any case could be closer to reality than thin matter sheets. 

Bonjour, Charmousis and Gregory have used the (2, 1; 1) version of the general split ()15() . (jlOj) . ()17() with 
scalar matter and consider both thick walls and the thin-wall limit perturbatively in the presence of gravity. 
They further specialize to the case of a spherical domain wall, which they show collapses. It would be worthwhile 
if this sort of example, which combines the modelling assumptions II) and III), is investigated using the 
(n, 0; -1) split subject to flat and to AdS asymptotics. 

For the particular example mentioned above, the matter profile is sigmoidal. However, to approximate a thin 
brane, one would want instead a hump profile, as indicated in Fig 9a. One could still use an inner b.c at z = 1, 
where now Kij = so that the linear Robin b.c 



dw (p 



(100) 



holds. One could in fact consider 2 problems (see Fig 9b): the flat sheet with directional asymptotics and the 
spherical sheet with a suitable inversion-in-the-sphere in place of the above inner b.c. Both of these problems 
avoid the main obstacles to existence and uniqueness proofs by possessing linear b.c's. The second problem is 
also invertible to an inner problem which is on a bounded region with no corners, for which results in [HEI apply- 

These problems are for 'plain' thick branes as opposed to branes containing spherically-symmetric objects. 
The first problem ought to more easily admit this extension. One might wish to study 'plain' thin branes too, 
e.g in the context of branes in relative motion. 

Finally, suppose instead one attempted to use Magaard's method to construct (4, 0) thin-matter-sheet data. 
Then one benefits from s = and the identification of xi = const with the z = 1 brane prescribes the topology, 
removes some of the sources of nonuniqueness and the brane's energy-momentum endows physical significance 
upon the (3, 0) covariance of the elimination procedure. In the thin matter sheet case, one is blocked by the 
nonanalyticity of p, but at least some thick matter sheet models could be built in this way. 



7 Discussion: Modelling Assumptions and Stringy Features 

We first address Q2 by discussing how the work in this paper may be generalized and which of these generaliza- 
tions may be regarded as stringy. We end with the conclusions of this paper as regards Ql . 

We did not consider the possibility of a higher codimension IHUI which could substantially alter the nature 
of the embeddings used. For models with 2 or more times, it is simply not possible to stick entirely to our 
suggestion to build on established mathematical physics. Finally, our study could be complicated by having 
more than the Einstein tensor in the bulk. For example, one can have a nontrivial Lovelock tensor term |91j in 
dimension 5, which arises from the variation of a Gauss-Bonnet term in the action, or one could have genuine 
higher-derivative terms [55]. Whilst our incipient results such as the correspondence between (n, 0; -1) and 
(n - 1, 1; 1) schemes at the simplest level or our BEFE ambiguity will have counterparts in these theories, these 
theories' equations are much larger than the EFE's, and their CP/IVP is much less well-studied than that of 
GR. Higher derivatives or non-standard bulk matter might lead to the violation of the energy conditions used in 
this paper. 

We next discuss which features of GR-based braneworld models are desirable for string theory. First for some 
caveats. We introduced the SMS braneworld as a GR generalization of the Randall-Sundrum scenario, moreover 
one that was as yet not sufficiently general for the purposes of GR in which the nature and existence of thin 
matter sheets at all times might not be a representative presupposition. It could be that it is other aspects of the 



Randall-Sundrum scenario that are of interest in string theory since it also has particle physics aspects to it and 
because it is a toy model of Hof ava-Witten theory. One may prefer to consider only such stringy scenarios rather 
than Randall-Sundrum or SMS. Also a great source of tension between string theorists and relativists is that the 
latter believe GR is suggestive of the importance of background-independent theories, whilst there is as yet no 
background-independent formulation of string theory available. Such ideas as conserved charges, flat sheets and 
simple bulks arise from assuming Minkowski or AdS backgrounds. Since GR is background-independent, one is 
interested there in generic solutions to the full field equations. This makes the notion of asymptotics necessary 
and subtle, and not necessarily compatible with flat sheets nor conserved charges. From a general perspective it 
is doubtful whether highly-specialized solutions such as Minkowski or AdS are likely to reflect the gravitational 
physics of the universe, unless it can be shown that large classes of generic solutions behave likewise or are 
attracted toward such solutions. 

String theory might favour many extra dimensions, although it is not clear how many of these should be large 
and how many compactified. Among the corrections to the Einstein tensor mentioned above, string theorists 
might favour the presence of the Gauss-Bonnet correction, since this is the only first-order correction for a 
heterotic string [HHI; an infinite series of correction terms is predicted. String theorists might favour very fiat 
walls on the grounds that these could be stabilized by being the lowest-energy carriers of BPS charge 0^5- These 
arguments would rule against the use of highly symmetric bulk manifolds/orbifolds, although string theorists 
might try to invoke them as fixed backgrounds so as to perturb about them. There may be different reasons for 
stability in GR and string theory. For example, in contrast to string-theoretic arguments favouring reflection- 
symmetric orbifolds, a GR domain wall could be stable precisely because the (3, 1) bulk is different on each 
side, as the wall separates domains in which symmetry is differently broken. On the whole, the string theorist 
would agree that bulks could be complicated by the presence of bulk matter from the closed string multiplet, as 
mentioned in Sec 4.5. At least in some models, gauge fields would be expected to occur only as fields confined to 
branes since they belong to the open string multiplet and some branes are where open strings end l3]. As strictly, 
the branes are of Planckian thickness, so it is not a disaster that thick branes are favourable toward rigorous 
mathematics for GR-based scenarios (this is recycling the argument in footnote 11). Of course, one could worry 
that they are thin enough for quantum gravitational effects to be important. 

We now address Ql. We feel that it is secure to say that GR-based braneworld models only make sense 
if realistic examples of them are stable to the dynamics which governs GR. We have investigated two schemes 
within which this study might be furthered: the (3, 1; 1) and (4, 0; -1) schemes. 

As regards (3, 1; 1) schemes, before pressing on with something like the SMS braneworld by perturbation 
theory and dynamical systems analysis, we think that it is important to check that the SMS braneworld's 
interpretation is properly understood. We have demonstrated that there exist many reformulations of SMS's 
BEFE's, thus making it clear that it is a choice as to whether the BEFE's contain SMS's Weyl term (requiring 
closure as some complicated third-order system such as in SMS's paper), and that it is a separate choice as to 
whether any terms homogeneously quadratic in the thin matter sheet energy-momentum tensor feature in the 
BEFE's. We have used this to illustrate the danger in truncating the BEFE's by setting bulk terms to be zero 
or to take a particular form. The full (4, l)-dimensional picture of the brane is required. Building this picture 
by sideways evolution appears unlikely to be well-behaved mathematically or to make sense from a physical 
(causal) point of view. We hope that this paper will either stimulate new research where these concepts are 
carefully thought about or force an admission that it is preferable to use a (4, 0; -1) scheme. We have argued 
that the removal of singularities by the (3, 1; 1) scheme is unlikely to be generally possible in a meaningful 
way. Whichever of the two schemes one works in, one is well-advised to heed the old, careful literature about 
solving the Gauss-Codazzi system rather than starting afresh. We have carefully detailed which parts of this 
traditionally (n, 0) scheme easily pass over to the (r, 1) case. 

As regards the (4, 0; -1) scheme that we have argued is favourable, we have argued that its IVP step is 
cleaner than its CP step and by itself serves to answer a number of questions. We have provided a survey of 
obstacles and ideas to extend this IVP with a thin matter sheet away from the relatively trivial time-symmetric 
cases hitherto considered. We also suggest that thick branes are more easily treated rigorously and that this 
choice should not affect the underlying physics if it is respectably robust. We feel that the next stage should be 
to solve for such data numerically. As this field grows, the underlying analytic mathematics will deserve detailed 
attention. 

We note that the scope of the (4, 0; -1) scheme would go as far as the study of colliding black holes, which 
are expected to be an important source of gravity waves. One requires the kind of methods that we promote in 
this paper for GR-based braneworlds for the full study of such high-curvature processes. Whilst this study is a 
tall order even for the usual GR j531 l2Uj . these calculations may be required if braneworld effects can be large 



enough to substantially alter the extreme physics of colliding black holes. 

Finally, we summarize the aspects of Q2 that should be asked as the partial answers to Ql begin to be 
built up. Does the introduction of stringy features change the outcome to any of issues? And so, are GR-based 
braneworlds adequate or typical frameworks for string theory? What scenarios should one explore in order to 
furnish string theory with reasonably unique predictions? 
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Figure captions 



Figure 1 

a) For unspecified or general signature, we use the following notation. Let (Af, a (rt + l)-dimensional spacetimc. Let dq be an embedded 
spacelike hypersurface at parameter value Aq, with unit normal n"^. The induced metric on E is hab = 9ab ~ Consider another nearby 

spacelike hypersurface Ti at parameter value Ai. Then nearby slices fit together in accordance with the figure, where the lapse is a = . ^ 

V Is I 

and the shift is j3i = g^i. In the particular case of a (n, 0;-l) split, we denote T by S, by and A by t whilst for a (n - 1, 1; 1) split we 
denote wc denote T by El, by and A by z. 

b) When we simultaneously talk of several splits in this paper, we use the notation presented in this figure. The white arrows represent the 
extrinsic curvatures associated with each of the cmbcddings. 

Figure 2 

A general means of study of hyperbolic p.d.e's is to require that the data be prescribed on spacelike surfaces [{n, 0; —1) scheme] a), not on 
timelike surfaces [(r, 1, 1) scheme], b). The line J in a) represents the inclusion of a snapshot of a thin matter sheet, whilst the whole of the 
surface in b) could be the worldsheet of an arbitrarily thin matter sheet. 

Figure 3 (Domain of dependence property) We should not claim to know too great a portion of the future, a) Given data for a hyperbolic 
system on a piece of a spacelike surface, we can predict the future only in the forward wedge (domain of dependence) within which all causal 
curves, i.e allowed paths of information fiow, emanate from the data. Outside this wedge, for all we know, large disturbances (represented in the 
picture by gravity waves) could influence the future arbitrarily soon. 

b) If we claim to know the 4-dimensional DOD, we preclude the influence by the extra dimension, thus making it largely redundant. The 
5-dimensional DOD of a thin sheet is small because of the thinness of the sheet. Thus gravity waves could be lurking very close in the bulk. This 
corresponds to the 4-d branewold EFE system not being closed. 
Figure 4 

a) The bucket-shaped construction that motivates the use of Sobolev spaces for the Cauchy problem for the Klein— Gordon equation. 

b) The analogous construction for the sideways Gauchy problem, which we argue is faulty on several counts. Sobolev spaces do not appear to be 

appropriate for the sideways Cauchy problem . 
Figure 5 

a) n is a (r, 1) (i.e timclikc) hypersurface containing a data strip T built by the Magaard method from prescribed 'data for the data' on a 
timclikc xi = 0. The shaded region Vl is the region of the (r + 1, 1) spacetimc supposedly entirely controlled by T i.e for which the Campbell 
method can validly provide the evolution. Typically T will not cover all of 11 in the xi-direction, so there will be nearby points like p outside T. 
Consider the future light-cone with apex p. No matter how thin in the z-diicction one considers £7 to be, it is always pierced by causal curves 
such as 7 on or in the light-cone, by taking 7 to be at a sufficiently slender angle (j> to Ft. Thus information can leak into f2 from elsewhere than 
P, which is a contradiction. Therefore, in parts at least, Cl is arbitrarily thin in the z-direction. One can envisage that sometimes there will be 
enough points like p that there is no region Q, at all. 

b) Suppose = is spacelike and is a section of the entire hypersurface n and T' is the data strip built from it by the Magaard method. There 
is now no room for points like p\ However, typically T' will not cover all of 11 in the a;i-direction, so there will be nearby points like q in the 
causal past J~{P') not in P'. No matter how thin one considers f2 to be, it is pierced by causal curves 7' from q which do not pass through P' 
, by taking these to be at a sufficiently slender angle cp' . So again, in parts at least, Q is arbitrarily thin and sometimes there will be enough 
points like q that there is no region Q at all. Small pieces of timelike hypersurfaces need not hold useful information for hyperbolic-type systems! 

Figure 6 

a) A simple illustration of how 4-dimensional expansion may be interpreted as being mainly due to shear in a 5-dimensional embedding manifold. 

b) Diagram of the embedding of flat FLRW universes into Minkowski spacetime in the standard 5-d Minkowski coordinates (T, Xi , X2 , X3 ,Z). 
The curved surfaces are the FLRW spacetimes. As one approaches T = (corresponding to t = in FLRW coordinates), each of these surfaces 

becomes tangent to the light cone (characteristics of Minkowski). The foliation by these surfaces also becomes bad here because the surfaces 
intersect. Note also that the FLRW geodesies move within each of the surfaces whereas the Minkowski geodesies clearly pierce these surfaces. 
Thus the 4-d and 5-d geodesies in this example are not the same. 
Figure 7 

The spread of the region of a snapshot P| "^Pl containing the junction J between the 2 spacelike bulk pieces E/^ and E^ upon 
which the prescribed data is smooth or analytic. Then the evolutions in the domains of dependence 15+ (E^,) and 15+ (E^) are smooth or analytic, 
but the causal future of J will typically be H^. 
Figure 8 

a) Excision of the inner region in (3, 0) black hole data construction by an inversion-in-the-sphere isometry identification between two copies of 
the black hole. This is along the lines of the method of images in electrostatics. 

b) Excision of the inner region in (4, 0) braneworld black hole data construction would involve mailing a guess r = f{z) for the excision region 
and then performing inversion-in-the-sphere isometries pointwise to identify two copies of the black hole. 

c) One is then to numerically solve for the data. One can then find the apparent horizon. If this intersects r = f{z) (horizon 1) then the data 
may be bad in that the part of the shaded portion of the excised region may be in causal contact with its surroundings. If the apparent horizon 
extends too far past the end of the excision region (horizon 2) , one might worry that the singularity also extends past the excised region in which 
case some of the data prescribed along r = should not have been prescribed since it lies on the singularity. 

Figure 9 

a) The metric, extrinsic curvature and matter profiles for thin branes and thick branes. 

b) Some of the proposed thick matter sheet IVP's; planar thick branes and spherical thick branes, which permit the use of a standard rather than 
directional notion of asymptotics. 



